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Abstract 



One may represent polynomials not only by their coefficients but also 
by arithmetic circuits which evaluate them. This idea allowed in the last 
fifteen years considerable complexity progress in effective polynomial equation 
solving. We present a circuit based computation model which captures all 
known symbolic elimination algorithms in effective Algebraic Geometry and 
exhibit a class of simple elimination problems which require exponential size 
circuits to be solved in this model. This implies that the known, circuit based 
elimination algorithms are already optimal. 
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1 Introduction 

We introduce and motivate a new computation model which is well-adapted to 
scientific computing in effective Algebraic Geometry and especially to elimination 
theory. This model is based on the symbolic manipulation of arithmetic circuits 
which evaluate rational functions. 
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problems. 



Most algorithms in effective Algebraic Geometry may be formulated as routines 
which operate on polynomials and rational functions and this suggests the represen- 
tation of these mathematical entities by arithmetic circuits. Thus, we shall consider 
arithmetic circuits as objects (in the sense of object oriented programming) which 
become mapped into abstract data types consisting of polynomials and rational 
functions. Inputs and outputs of the routines of our computation model will there- 
fore be arithmetic circuits. These circuits will appear as parameterized in the sense 
that they depend on two distinct ingredients, namely on piecewise rational func- 
tions called "basic parameters" and indeterminates, called "input variables". On 
the other hand, so called "elementary routines" will constitute the basic building 
block of our computation model. They will be branching-free and therefore our 
parameterized arithmetic circuits will be branching-free too. 

However, in effective Algebraic Geometry, divisions are sometimes unavoidable 
and divisions may lead to branchings. Nevertheless, in typical situations, they may 
be replaced by limit processes. In order to capture this situation, we introduce in 
Section |3~T the notion of a robust parameterized arithmetic circuit. 

An important issue will be the concept of well behavedness of routines, under 
certain modifications of the input circuits. This concept will appear in Section 
|3.3| in different disguises, called well behavedness under restrictions and reductions, 
(output) isoparametricity and coalescence. 

All these technical notions have in common that they allow to formulate algo- 
rithmical restrictions on the routines of our computation model which are motivated 
by specific quality attributes of programs in Software Engineering. 

In this sense we establish first in Section 3.3.2 a branching-free variant of our 
computation model which operates on robust parameterized arithmetic circuits. The 
algorithms captured by this variant are the elementary routines mentioned before. 

In order to capture the whole spectrum of really existing elimination algorithms 
in Algebraic Geometry, we extend our computation model in Section 3.3.3 admit- 
ting some limited branchings. The resulting algorithms are called branching par- 
simonious. Moreover we introduce the concept of a procedure as a branching par- 
simonious algorithm with a particular architecture. Procedures are well suited to 
discuss computational issues in effective elimination theory. In Section we apply 
our computation model to this task. 

Before we are going to enter into the details of these applications we look closer 
to our computation model. The basic construction method of elementary routines is 
recursion. Mimicking the directed acyclic graph structure of a given robust param- 
eterized arithmetic circuit we compose the graphs of previously fixed computations 
in order to obtain a new parameterized arithmetic circuit. To guarantee that the 
resulting circuit is again robust, we require that these compositions behave well un- 
der restrictions. Moreover, we require that at corresponding nodes the intermediate 
results of the two circuits become directly linked by a continuous, piecewise rational 
map and call this requirement isoparametricity. This means that our computation 



model includes also formal specifications. We develop this aspect in Sections 3.1.2 
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and 13.3.2.41 

In Section [4] we use our computation model to show that already very elementary 



elimination problems require exponential time for their solution (see Theorem 10 



Proposition 11 and Theorem 12 below). 



In particular, we exhibit in Section |4T3| a family of parameterized boolean circuits 
whose (standard) arithmetizations represent an elimination problem which requires 
exponential time to be solved in our model (see Theorem 13 below). 

As a major outcome of this paper we exhibit in Section 4.5 an infinite family 
of parameter dependent elimination polynomials which require essentially division- 
free, robust parameterized arithmetic circuits of exponential size for their evaluation, 
whereas the circuit size of the corresponding input problems grows only polynomi- 
ally. We observe that essentially division-free, robust parameterized arithmetic 
circuits for elimination polynomials capture the intuitive meaning of an algorithmic 
solution with few equations and branchings of the underlying elimination problem. 

The proof of this result, which is absolutely new in his kind, is astonishly ele- 
mentary and simple. 

In Section 14.61 we arrive at the conclusion that our method to show lower com- 
plexity bounds consists of counting how many steps are necessary to decompose a 
given rational map into a sequence of "simple" blow ups and a polynomial map. 

Finally in Section 4.7 we establish a link between our computation model and our 
lower bound results with other complexity views in geometric elimination theory. 
In this context we discuss the BSS-model of [BSS89| and the view of interactive 
protocols. 

Our computation model and complexity results are based on the concept of a 
geometrically robust constructible map. This concept was introduced in [GHMSllJ 
and we develop it further in Section [2j which is devoted to the algebraic geometric 
underpinning of the present paper. 

The relevance of the lower complexity bounds of this paper for elimination prob- 
lems depends on the "naturalness" of the computation model. Therefore we empha- 
size throughout this article the arguments which justify our computation model. Of 
course, these arguments cannot be entirely of mathematical nature. In this paper 
they are borrowed from Software Engineering which constitutes a discipline which 
analyzes and qualifies practical programming issues. In these terms we show that a 
circuit based algorithm which solves most elementary parametric elimination prob- 
lems and which is programmed under the application of the most common rules of 
Software Engineering, can never be efficient. 

This paper is based on the idea to represent, e.g., in the case of elimination 
algorithms, polynomials by arithmetic circuits, which are not considered as algo- 
rithms but as data structures. This has sometimes certain advantages. For ex- 
ample, the generic n x m determinant has as polynomial in the matrix entries n\ 
terms but it can be evaluated by division and branching-free arithmetic circuit of 
size 0(n 5 ) (|Mul87]). More generally, an arbitrary elimination polynomial of degree 
S and m variables has always a circuit representation which is essentially of order 
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8°^\ whereas the representation by its coefficients may become of order S n ^ (see 
|GHM+98j and |GHH+97| for details). 

The idea of the circuit representation of polynomials was introduced in The- 
oretical Computer Science at the beginning of the eighties by the first author in 
collaboration with Malte Sieveking (Frankfurt/Main). However, the first publica- 
tions on this subject treated only the case of the elimination of a single variable (see 
e.g. [HS8T] . |Kal88j and |FIK86] ). 

In the case of the simultaneous elimination of several variables, substantial 
progress was made in the nineties by the first author in collaboration with Marc 
Giusti (Paris) and Jacques Morgenstern (Nice). A series of paper give account of 
this development (see e.g. |GH93| . |HM93| . [GHM+98] and |GHH+97j ). In partic- 
ular the work with Jacques Morgenstern lead to new views which influenced this 
paper. 



2 Concepts and tools from Algebraic Geometry 



In this section, we use freely standard notions and notations from Commutative Al- 
gebra and Algebraic Geometry. These can be found for example in |Lan93j . [ZS60J, 



|Kun85j and |Sha94] . In Sections 2.2 and 2.3, we introduce the notions and defini- 
tions which constitute our fundamental tool for the modelling of elimination prob- 
lems and algorithms. Most of these notions and their definitions are taken from 
jGHMSllj . 



2.1 Basic notions and notations 

For any n e N, we denote by A™ := A n (C) the n-dimensional affine space C n 
equipped with its respective Zariski and Euclidean topologies over C. In algebraic 
geometry, the Euclidean topology of A n is also called the strong topology. We shall 
use this terminology only exceptionally. 

Let X±, . . . , X n be indeterminates over C and let X := (Xi, . . . , X n ). We denote 
by C[X] the ring of polynomials in the variables X with complex coefficients. 

Let V be a closed affine subvariety of A n . As usual, we write dim^ for the 
dimension of the variety V. Let C\, . . . , C s be the irreducible components of V. For 
1 < J < s we define the degree of Cj as the number of points which arise when 
we intersect Cj with dimCj many generic affine hyperplanes of A n . Observe that 
this number is a well-determined positive integer which we denote by deg Cj. The 
(geometric) degree degV of V is defined by degV := Xa<j< s degCj. This notion 
of degree satisfies the so called Bezout Inequality. Namely, for another closed affine 
subvariety W of A n we have deg V fl W < deg V ■ deg W. 

For details we refer to [Hci83j, where the notion of geometric degree was intro- 
duced and the Bezout Inequality was proved for the first time (other references are 
jFul84] and |Vog84| ). 
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For fx, . . . , f s , g G C[X] we shall use the notation {fi = 0, . . . , f s = 0} in order 
to denote the closed affine subvariety V of A n defined by fi, . . . , f s and the notation 
{/i = 0, . . . , f s = 0,g 7^ 0} in order to denote the Zariski open subset V g of V 
defined by the intersection of V with the complement of {g = 0}. Observe that V g 
is a locally closed affine subvariety of A™ whose coordinate ring is the localization 
C[V] g oiC[V}. 

We denote by I(V) : = {/ G C[X] : f(x) = for any x G V} the ideal of definition 
of V in C[X] and by C[V] := : V ->■ C ; there exists / G C[X] with p(z) = 
/(x) for any x G its coordinate ring. Observe that C[V] is isomorphic to the 
quotient C-algebra C[V] := C[X]/I(V). If V is irreducible, then C[V] is zero- 
divisor free and we denote by C(V) the field formed by the rational functions of V 
with maximal domain (C(V) is called the rational function field of V). Observe that 
C(V) is isomorphic to the fraction field of the integral domain C[V]. 

In the general situation where V is an arbitrary closed affine subvariety of A n , 
the notion of a rational function of V has also a precise meaning. The only point to 
underline is that the domain, say U, of a rational function of V has to be a maximal 
Zariski open and dense subset of V to which the given rational function can be 
extended. In particular, U has a nonempty intersection with any of the irreducible 
components of V. 

As in the case where V is irreducible, we denote by C(V) the C-algebra formed 
by the rational functions of V. In algebraic terms, C(V) is the total quotient ring 
of C[V] and is isomorphic to the direct product of the rational function fields of the 
irreducible components of V. 

Let be given a partial map : V --- » W, where V and W are closed sub varieties 
of some affine spaces A n and A m , and let ^i, ... , <fi m be the components of <fi. With 
these notations we have the following definitions: 

Definition 1 (Polynomial map) The map <fi is called a morphism of affine va- 
rieties or just polynomial map if the complex valued functions <pi, . . . , <p m belong to 
C[V}. Thus, in particular, (ft is a total map. 

Definition 2 (Rational map) We call <p a rational map ofV to W , if the domain 
U of <f) is a Zariski open and dense subset of V and 0i, . . . , <p m are the restrictions 
of suitable rational functions of V to U. 

Observe that our definition of a rational map differs from the usual one in Al- 
gebraic Geometry, since we do not require that the domain U of <f> is maximal. 
Hence, in the case m := 1, our concepts of rational function and rational map do 
not coincide (see also |GHMSll] ). 

2.1.1 Constructible sets and constructible maps 

Let A4 be a subset of some affine space A™ and, for a given nonnegative integer m, 
let : M. — ■* A m be a partial map. 
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Definition 3 (Constructible set) We call the set Ai constructible if Ai is de- 
finable by a Boolean combination of polynomial equations. 

A basic fact we shall use in the sequel is that if Ai is constructible, then its 
Zariski closure is equal to its Euclidean closure (see, e.g., |Mum88j . Chapter I, §10, 
Corollary 1). In the same vein we have the following definition. 

Definition 4 (Constructible map) We call the partial map constructible if the 

graph of (p is constructible as a subset of the affine space A n x A' m . 

We say that is polynomial if is the restriction of a morphism of affine varieties 
A n A m to a constructible subset Ai of A n and hence a total map from Ai to A m . 
Furthermore, we call a rational map of Ai if the domain U of is contained in 
Ai and is the restriction to Ai of a rational map of the Zariski closure Ai of Ai. 
In this case U is a Zariski open and dense subset of Ai. 

Since the elementary, i.e., first-order theory of algebraically closed fields with 
constants in C admits quantifier elimination, constructibility means just elementary 
definability. In particular, is constructible implies that the domain and the image 
of are constructible subsets of A™ and A" 1 , respectively. 

Remark 1 A partial map : Ai —* A m is constructible if and only if it is piecewise 
rational. If is constructible there exists a Zariski open and dense subset U of Ai 
such that the restriction (f)\u of to U is a rational map of Ai (and of Ai ). 

For details we refer to [GHMSll], Lemma 1. 

2.2 Weakly continuous, strongly continuous, topologically 
robust and hereditary maps 

We are now going to present the notions of a weakly continuous, a strongly con- 
tinuous, a topologically robust, a geometrically robust and a hereditary map of a 
constructible set Ai. These five notions will constitute our fundamental tool for the 
modelling of elimination problems and algorithms. 

Definition 5 Let Ai be a constructible subset of A n and let : Ai — > A m be a 
(total) constructible map. We consider the following four conditions: 

(i) there exists a Zariski open and dense subset U of Ai such that the restriction 
<j)\u of to U is a rational map of Ai and the graph of is contained in the 
Zariski closure of the graph of <p\u in Ai x A m ; 

(ii) is continuous with respect to the Euclidean, i.e., strong, topologies of Ai and 
A m ; 
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(Hi) for any sequence (xk)km of points of M. which converges in the Euclidean 
topology to a point of M., the sequence (4>(xk))keN is bounded; 

(iv) for any constructible subset M of Ai the restriction : M — > A m is an 
extension of a rational map of M and the graph of <f>\j^ is contained in the 
Zariski closure of the graph of this rational map in M x A m . 

We call the map <f> 

• weakly continuous if <p satisfies condition (i), 

• strongly continuous if (p satisfies condition (ii), 

• topologically robust if <ft satisfies conditions (i) and (Hi), 

• hereditary if <fi satisfies condition (iv). 

In all these cases we shall refer to Ai as the domain of definition of or we shall 
say that (ft is defined on Ai. 

Lemma 2 ( jGHMSll^ . Lemma 4) A strongly continuous constructible map is al- 
ways weakly continuous, topologically robust and hereditary. 



In Section 2.3, we shall establish an algebraic condition, namely geometric ro- 



bustness, which implies hereditarity. 

2.3 The concept of robustness for constructible maps 

In this Section we introduce the algebraic-geometric tools we shall use in Section 
[3] and [4] for the mathematical modelling of algorithms which solve parameterized 
computational problems. The main issue of this section will be the notion of a 
geometrically robust constructible map which captures simultaneously the concepts 



of topological robustness and hereditarity introduced in Section |2.2 

We first characterize in algebraic terms the concept of topological robustness 
(Theorem [3] below). In Section [3] we shall interpret topological robustness as the 
informal concept of coalescence (we call it informal because distinct authors intro- 
duce it differently, following the context). For example in Interpolation Theory 
coalescence refers to certain types of "convergence" of problems and algorithms (see 
|BC97j . |dBR92| . |Olv06] and |GHMS11| for details). In this paper coalescence will 
be the algorithmic counterpart of topological robustness. 

Finally, we introduce the notion of a geometrically robust constructible map and 
show that such maps are always hereditary. In particular they are topologically 
robust and give rise to coalescent algorithms. 
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2.3.1 An algebraic characterization of the notion of topological 
robustness 



In this subsection, we present an algebraic-geometric result of [OHMSllJ which will 
be relevant in Sections 2.3.2[ [3] and [|| 



For the moment let us fix a constructible subset Ai of the affine space A n and a 
(total) constructible map (p : Ai — > A' m with components (pi, . . . , (p m . 

We consider now the Zariski closure Ai of the constructible subset Ai of A n . 
Observe that Ai is a closed affine subvariety of A n and that we may interpret C(Ai) 
as a C[.M] -module (or C[.M] -algebra). 

Fix now an arbitrary point x of Ai. By 9Jl x we denote the maximal ideal of 
coordinate functions of C[«M] which vanish at the point x. By CLVl]^ we denote 
the local C-algebra of the variety Ai at the point x, i.e., the localization of C[Af] at 
the maximal ideal Wl x . By C(Ai)<m x we denote the localization of the C[.M] -module 

c(M) at m x . 

Following Remark [TJ we may interpret (pi, ... , (p m as rational functions of the 
affine variety Ai and therefore as elements of the total fraction ring C(Ai) of C[jVf]. 
Thus C[yVf][0i, . . . , 4>m] and C[Ai]gjr a [^ij • • • > 4>m\ are C-subalgebras of C(Ai) and 
C(Ai)yi x which contain C[.M] and CfA^]^, respectively. 

With these notations we are able to formulate the following statement which 
establishes the bridge to an algebraic understanding of the notion of topological 
robustness. 

Theorem 3 (' L GHMSlT$, Corollary 11) Let notations and assumptions be as be- 
fore and suppose that the constructible map : Ai — > A m is weakly continuous. 
Then <fi is topologically robust if and only if for any point x of At the C-algebra 
C[At}<xn x [4>i, . . . , 4>m\ is a finite C[At]<xn x -module. 

The only if part of Theorem [3] is an almost immediate consequence of |CGH + 03j , 
Lemma 3, which in its turn is based on a non-elementary and deep result from 
Algebraic Geometry, namely Zariski's Main Theorem (see, e.g., |Ive73] . §IV.2). 

Let 4> : Ai — > A m be a topologically robust constructible map and let u be an 
arbitrary point of Ai. From Theorem [3] one deduces easily that for all sequences 
(uk)keN of points Uk € Ai which converge to u, the sequences (0(«fc))fceN have only 
finitely many accumulation points. 



2.3.2 The notion of geometrical robustness 

The main mathematical tool of Section [3] of this paper is the notion of geometrical 
robustness we are going to introduce now. We shall use the same notations as in 
Section 12.3.11 

Definition 6 Let Ai be a constructible subset of a suitable affine space and let 
<p : Ai — > A m be a (total) constructible map with components (pi, ... , (p m . According 
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to Remark [i| we may interpret 0i, . . . , m as rational functions of Ai. We call cf) 
geometrically robust if for any point the following two conditions are satisfied: 

(i) ClyW]^ • • • ) 0m] is a finite C[Ai]<m x -module. 

(ii) C[Ai]<xn x [(f>i, . . . , <p m ] is a local C[Ai]<xn x -algebra whose maximal ideal is gener- 
ated by ffl x and (fix - <pi(x), . . . ,<f) m - <fi m (x). 

Observe that the notion of a geometrically robust map makes also sense when C 
is replaced by an arbitrary algebraically closed field (of any characteristic). In view 
of Theorem [3] the same is true for the notion of a topologically robust map. 

In this paper we shall restrict our attention to the algebraically closed field C. In 
this particular case we have the following characterization of geometrically robust 
constructible maps. 

Theorem 4 Let assumptions and notations be as before. Then the constructible 
map : .M — 7- A m is geometrically robust if and only if <ft is strongly continuous. 

Proof. Suppose that the constructible map is geometrically robust. We are 
first going to show that <fi is weakly continuous. 

By Remark [T] there exists a Zariski open and dense subset U of Ai such that 
the restriction map <p\u is rational. Let Yi,...,Y m be new indeterminates, Y := 
(Yi, . . . ,Y m ) and suppose that the affine ambient space of Ai has dimension n. 
Observe that any (n + m)-variate polynomial over C which vanishes on the graph of 
the rational map <f>\u gives rise to a polynomial A G C[AT] [Y] with A[(f>i, . . . , m ] = 0. 

Let x be an arbitrary point of Ai and consider A as an element of CfA^]^ [Y]. 
Denote by A(x,<fi(x)) the value of A at (x,<fi(x)). Then condition (ii) of Defi- 
nition [6] implies that A[<fix, . . . , m ] — A(x, 0(x)) belongs to the maximal ideal of 
Ct-M]^^!, . . . , m ]. From A[<f>i, ... , m ] = we deduce now A(x, (j)(x)) = 0. 

Since the choice of x G Ai was arbitrary, we conclude that A vanishes on the 
graph of <p. This implies that the graph of is contained in the Zariski closure of 
the graph of (f)\jj. Hence <fi is weakly continuous. 

Let be given an arbitrary point x 6 Ai and a sequence (xk)km, x k G Ai, which 
converges to x in the strong topology of Ai. We are now going to show that the 
sequence ((f)(xk))ken converges to 4>(x). 

Since <fi is weakly continuous, we deduce from condition (i) of Definition [6] and 
Theorem [3] that the sequence (0(x n ))fceN contains at least one accumulation point, 
say a = (ax, ■ ■ ■ , a m ), which belongs to A m . Let a be the ideal of all polynomials 
A G CLVl]^ [Y] that vanish at the point (x, a) G A n x A m . Without loss of generality 
we may assume that the sequence ((p(xk))keN converges to a. Let a := {^4(0); A G a} 
be the image of the ideal a under the surjective C[A / f]<H a; _a lg e t ) i'a homomorphism 
C[M]m x [Y]^ -> C[Ai] mx [0i, LLJ _, m ] which maps Yi, . . . , Y m onto 0i, . . . , m . Ob- 
serve that a is an ideal of CfA^]^ [0i, . . . , <p m ). 

We are now going to show the following statement. 

Claim 5 The ideal a is proper. 
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Proof of the claim. Suppose that the ideal a is not proper. Then there exists 

a polynomial A = £\ ", -,- )y: ■ ■ of a > with a h-j m e C[M]m x , which 

satisfies the condition j m a h-j m 4^i ■ ■ ■ 0m = A(<fi) = 1. Since for any m- 

tuple of indices ji, . . . , j m the rational function aj 1 ,,,j m of .M is defined at x and 
the sequence (xk)k<m converges to x, we may assume without loss of generality that 
aj x _j m is defined at x^ for any k G N and that (%i...j m O E fc))fceN converges to cij 1 ...j m {x). 
We may therefore write A^ x > := a h...j m { x ')Yi 1 ■ ■ ■ ^m m ^ C[F] for x' := x or 
x' := Xfc, k G N. From A G a we deduce A^'(a) = 0. By assumption (0(xfc))k 6 N 
converges to a. Hence the sequence of complex numbers (A^ Xk ^ ((j)(xk))) k&n converges 
to A^ x '(a) = 0. On the other hand A((j>) = 1 and the weak continuity of <p imply 
A( Xk \(j)(xk)) = 1 for any k G N. This contradiction proves our claim. 

From condition (ii) of Definition [6] we deduce that the CfA'ljgjt^^algebra 
C[A^]ot x [01, • • • > 0m] contains a single maximal ideal, say Wl, and that VJt is generated 
by Wl x and 0x - 0i(x), . . . , m -J> m {x). 

Since by Claim [5] the ideal a is proper, a must be contained in Observe 
that the polynomials Y\ — a±, . . . , Y m — a m belong to a. Hence </>i — ai, . . . , 4> m — d m 
belong to a and therefore also to VJt. Since 071 is proper, this is only possible if 
ai = 4>i{x), . . . , a m = 4>m{x) holds. 

Thus the sequence (0(xfc))fe£N converges to 0(x). 

Suppose now that the constructible map is strongly continuous. From Lemma [2] 
we deduce that is topologically robust. Theorem [3] implies now that <fi satisfies 
condition (i) of Definition [6] at any point of M.. 

Let x be an arbitrary point of M.. We have to show that (j) satisfies at x condition 
(ii) of Definition [6j 

Since the graph of <ft is constructible, its strong and Zariski closures in M. x A m 
coincide. Moreover, since is by assumption strongly continuous, its graph is closed 
with respect to the strong topology of M. x A m and therefore also with respect to 
the Zariski topology. Let a be an arbitrary maximal ideal of the C [M] <xr x -algebra 
CLVfjajtaf^i, ■ ■ , 0m]- Then there exists a point a = (ai, . . . , a m ) of A m such that 
a is generated by Wl x and 0i — ai, . . . , TO — a m . Thus (x, a) G M. x A m belongs 
to the Zariski closure of the graph of in M. x A m and therefore to the graph 
of itself. This implies a = <p(x). With other words, a is generated by Wl x and 
0i — 0i (x), . . . , m — 4> m (x). There is exactly one ideal of CfA^]^ [<pi, . . . , ip m ] which 
satisfies this condition. Therefore the CfM]^ -algebra CfA^]^ [<fi, ■ ■ ■ ,(p m ] is local 
and condition (ii) and Definition [6] is satisfied at the point x G M. ■ 

Theorem [4] implies immediately the following result which will be fundamental 
in the sequel. 

Corollary 6 Geometrically robust constructible maps are weakly continuous, hered- 
itary and even topologically robust. If we restrict a geometrically robust constructible 
map to a constructible subset of its domain of definition we obtain again a geometri- 



10 



cally robust map. Moreover the composition and the cartesian product of two geomet- 
rically robust constructible maps are geometrically robust. The geometrically robust 
constructible functions form a commutative C-algebra which contains the polynomial 
functions. 

Notice that Corollary [6] remains mutatis mutandis true if the notions of geomet- 
rical and topological robustness are applied to constructible maps defined over an 
arbitrary algebraically closed field k. 

Adapting the corresponding proofs to this more general situation, one sees that 
weak continuity and hereditarity of geometrically robust constructible maps with 
irreducible domains of definition follows from [GHMSllJ, Proposition 16, Theorem 
17 and Corollary 18. These results imply also that restrictions of such maps to 
irreducible constructible subsets of their domains of definition are again geometri- 
cally robust. From this one deduces immediately the same statements for the case 
of arbitrary domains of definition. Topological robustness is a direct consequence 
of Definition |6j Closedness under composition follows from the transitivity law for 
integral dependence. One infers from Definition [6] closedness under cartesian prod- 
ucts and that the geometrically robust constructible functions form a commutative 
/c-algebra which contains the polynomial functions. 

Theorem [4] is new. It gives a topological motivation for the rather abstract, 
algebraic notion of geometrical robustness. The reader not acquainted with commu- 
tative algebra may just identify the concept of geometrical robustness with that of 
strong continuity of constructible maps. 

The origin of the concept of a geometrically robust map can be found, implic- 
itly, in [GH01J. It was introduced explicitly for constructible maps with irreducible 
domains of definition in |GHMSlT] , where it is used to analyze the complexity char- 
acter of multivariate Hermite-Lagrange interpolation. 

3 A software architecture based model for computations 
with parameterized arithmetic circuits 

3.1 Parameterized arithmetic circuits and their semantics 

The routines of our computation model, which will be introduced in Section |3.3 
operate with circuits representing parameter dependent rational functions. They will 
behave well under restrictions. In this spirit, the objects of our abstract data types 
will be parameter dependent multivariate rational functions over C, the concrete 
objects of our classes will be parameterized arithmetic circuits and our abstraction 
function will associate circuits with rational functions. In what follows, C may 
always be replaced, mutatis mutandis, by an arbitrary algebraically closed field (of 
any characteristic). 

Let us fix natural numbers n and r, indeterminates X%, . . . , X n and a non-empty 
constructible subset Ai of A r . By 7Ti, . . . , 7iy we denote the restrictions to Ai of the 
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canonical projections A r — >■ A . 

A (by Ai ) parameterized arithmetic circuit (with basic parameters 7Ti, . . . ,ir r 
and inputs X\, . . . , X n ) is a labelled directed acyclic graph (labelled DAG) satisfying 
the following conditions: 

each node of indegree zero is labelled by a scalar from C, a basic parameter 71*1, . . . , 7iy 
or a input variable X\, . . . , X n . Following the case, we shall refer to the scalar, basic 
parameter and (standard) input nodes of 0. All other nodes of have indegree 
two and are called internal. They are labelled by arithmetic operations (addition, 
subtraction, multiplication, division). A parameter node of depends only on scalar 
and basic parameter nodes, but not on any input node of (here "dependence" refers 
to the existence of a connecting path). An addition or multiplication node whose 
two ingoing edges depend on an input is called essential. The same terminology is 
applied to division nodes whose second argument depends on an input. Moreover, 
at least one circuit node becomes labelled as output. Without loss of generality we 
may suppose that all nodes of outdegree zero are outputs of 0. 

We consider as a syntactical object which we wish to equip with a certain 
semantics. In principle there exists a canonical evaluation procedure of assigning 
to each node a rational function of Ai x A n which, in case of a parameter node, may 
also be interpreted as a rational function of Ai. In either situation we call such a 
rational function an intermediate result of 0. 

The evaluation procedure may fail if we divide at some node an intermediate 
result by another one which vanishes on a Zariski dense subset of a whole irreducible 
component of Ai x A n . If this occurs, we call the labelled DAG inconsistent, 
otherwise consistent. From |CGH + 03j . Corollary 2 (compare also |HS80j . Theorem 
4.4 and [GH01J, Lemma 3) one deduces easily that testing whether an intermediate 
result of vanishes on a Zariski dense subset of a whole irreducible component of 
Ai x A n can efficiently be reduced to the same task for circuit represented rational 
functions of Ai (the procedure is of non-uniform deterministic or alternatively of 
uniform probabilistic nature). 

Mutatis mutandis the same is true for identity checking between intermediate 
results of 0. If Ai is irreducible, both tasks boil down to an identity-to-zero test 
on jM. In case that Ai is not Zariski dense in A r , this issue presents a major 
open problem in modern Theoretical Computer Science (see [Sax09j and |ShplO] for 
details). 

If nothing else is said, we shall from now on assume that is a consistent 
parameterized arithmetic circuit. The intermediate results associated with output 
nodes will be called final results of 0. 

We call an intermediate result associated with a parameter node a parameter of 
and interpret it generally as a rational function of Ai. A parameter associated 
with a node which has an outgoing edge into a node which depends on some input 
of is called essential. In the sequel we shall refer to the constructible set Ai as 
the parameter domain of 0. 
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We consider as a syntactic object which represents the final results of 0, i.e., 
the rational functions of Ai x A n assigned to its output nodes. In this way becomes 
introduced an abstraction function which associates with these rational functions. 
This abstraction function assigns therefore to a rational map Ai x A™ --- » A 9 , 
where g is the number of output nodes of 0. On its turn, this rational map may also 
be understood as a (by Ai) parameterized family of rational maps A n — > A q . 

Now we suppose that the parameterized arithmetic circuit has been equipped 
with an additional structure, linked to the semantics of 0. We assume that for 
each node p of there is given a total constructible map Ai x A™ — > A 1 which 
extends the intermediate result associated with p. In this way, if has K nodes, we 
obtain a total constructible map Q : Ai x A n — > A K which extends the rational map 
Ai x A™ --- > A K given by the labels at the indegree zero nodes and the intermediate 
results of 0. 

Definition 7 (Robust circuit) Let notations and assumptions be as before. The 
pair (0, Q) is called a robust parameterized arithmetic circuit if the constructible map 
Q is geometrically robust. 

We shall make the following two observations to this definition. 

We state our first observation. Suppose that (0, fl) is robust. This means that the 
constructible map fl : Ai x A™ — > A K is geometrically and hence also topologically 
robust and hereditary. Moreover, the above rational map M x A" — ■» A A can be 
extended to at most one geometrically robust constructible map Q : Ai x A n — > A K . 
Therefore we shall apply from now on the term "robust" also to the circuit 0. 

Let us now state our second observation. We may consider the parameterized 
circuit as a program which solves the problem to evaluate, for any sufficiently 
generic parameter instance the rational map A n --- » A q which we obtain by 

specializing to the point u the first argument of the rational map Ai x A n --- > A q 
defined by the final results of 0. In this sense, the "computational problem" solved 
by is given by the final results of 0. 

Being robust becomes now an architectural requirement for the circuit and for 
its output. Robust parameterized arithmetic circuits may be restricted as follows: 

Let Af be a constructible subset of Ai and suppose that (0, fl) is robust. Then 
Corollary [6] implies that the restriction £l\j\r X A n of the constructible map to Afx A n 
is still a geometrically robust constructible map. 

This implies that (0, fl) induces a by TV parameterized arithmetical circuit 0j^. 
Observe that 0^ may become inconsistent. If 0^ is consistent then (0m, ^IaTxA™) 
is robust. The nodes where the evaluation of 0^ fails correspond to divisions of 
zero by zero which may be replaced by so called approximative algorithms having 



unique limits (see Section 3.3.2). These limits are given by the map fl^ X A n - We 
call (0m, fl\jVxA n ), or simply 0n, the restriction of (0, Q) or to Af. 

We say that the parameterized arithmetic circuit is totally division-free if any 
division node of corresponds to a division by a non-zero complex scalar. 



13 



We call essentially division-free if only parameter nodes are labelled by divi- 
sions. Thus the property of being totally division-free implies that is essentially 
division-free, but not vice versa. Moreover, if is totally division-free, the rational 
map given by the intermediate results of is polynomial and therefore a geometri- 
cally robust constructible map. Thus, any by Ai parameterized, totally division-free 
circuit is in a natural way robust. 

In the sequel, we shall need the following elementary fact. 

Lemma 7 Let notations and assumptions be as before and suppose that the pa- 
rameterized arithmetic circuit is robust. Then all intermediate results of are 
polynomials in Xi, . . . ,X n over the C-algebra of geometrically robust constructible 
functions defined on A4 . 

Proof. Without loss of generality we may assume that A4 is irreducible. Let p be a 
node of which computes the intermediate result G p : M. x A n — > A 1 . Definition 6 
(z) and the irreducibility of M. imply that G p is a polynomial of C(A4)[X 1 , . . . , X n }. 
Observe that any x G A ra induces a geometrically robust constructible map M. — > A 1 
whose value at the point u G M. is G p (u,x). Using interpolation at suitable points 
of A™, we see that the coefficients of the polynomial G p are geometrically robust 
constructible functions with domain of definition M.. m 

The statement of this lemma should not lead to confusions with the notion of 
an essentially division-free parameterized circuit. We say just that the intermediate 
results of are polynomials in Xi, . . . , X n and do not restrict the type of arithmetic 
operations contained in (as we did defining the notion of an essentially division- 
free parameterized circuit). 

Whether a division of a polynomial by one of its factors may always be substi- 
tuted efficiently by additions and multiplications is an important issue in Theoretical 
Computer Science (compare |Str73j ). 

To our parameterized arithmetic circuit we may associate different complexity 
measures and models. In this paper we shall mainly be concerned with sequential 
computing time, measured by the size of 0. Here we refer with "size" to the number 
of internal nodes of which count for the given complexity measure. Our basic 
complexity measure is the non-scalar one (also called Ostrowski measure) over the 
ground field C. This means that we count, at unit costs, only essential multiplica- 
tions and divisions (involving basic parameters or input variables in both arguments 
in the case of a multiplication and in the second argument in the case of a division), 
whereas C-linear operations are free (see |BCS97j for details). 

3.1.1 Operations with robust parameterized arithmetic circuits 

3.1.1.1 The operation join 

Let 71 and 72 be two robust parameterized arithmetic circuits with parameter 
domain Ai and suppose that there is given a one-to-one correspondence A which 
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identifies the output nodes of 71 with the input nodes of 72 (thus they must have 
the same number). Using this identification we may now join the circuit 71 with the 
circuit 72 in order to obtain a new parameterized arithmetic circuit 72 *a 7i with 
parameter domain A4. The circuit 72 *a 7i has the same input nodes as 71 and the 
same output nodes as 72 and one deduces easily from Lemma [7] and Corollary [6] 
that the circuit 72 *a 7i is robust and represents a composition of the rational maps 
defined by 71 and 72, if 72 *a 7i is consistent. The (consistent) circuit 72 *a 7i is 
called the (consistent) join of 71 with 72. 

Observe that the final results of a given robust parameterized arithmetic circuit 
may constitute a vector of parameters. The join of such a circuit with another robust 
parameterized arithmetic circuit, say 0, is again a robust parameterized arithmetic 
circuit which is called an evaluation of (3. Hence, mutatis mutandis, the notion of 
join of two routines includes also the case of circuit evaluation. 

3.1.1.2 The operations reduction and broadcasting 

We describe now how, based on its semantics, a given parameterized arithmetic 
circuit (3 with parameter domain A4 may be rewritten as a new circuit over A4 
which computes the same final results as (3. 

The resulting two rewriting procedures, called reduction and broadcasting, will 
neither be unique, nor generally confluent. To help understanding, the reader may 
suppose that there is given an (efficient) algorithm which allows identity checking 
between intermediate results of (3. However, we shall not make explicit reference to 
this assumption. We are now going to explain the first rewriting procedure. 

Suppose that the parameterized arithmetic circuit /3 computes at two different 
nodes, say p and p', the same intermediate result. Assume first that p neither 
depends on p', nor p' on p. Then we may erase p' and its two ingoing edges (if p' is 
an internal node) and draw an outgoing edge from p to any other node of which is 
reached by an outgoing edge of p'. If p' is an output node, we label p also as output 
node. Observe that in this manner a possible indexing of the output nodes of (3 may 
become changed but not the final results of (3 themselves. 

Suppose now that p' depends on p. Since the DAG (3 is acyclic, p does not 
depend on p'. We may now proceed in the same way as before, erasing the node p'. 

Let (3' be the parameterized arithmetic circuit obtained, as described before, by 
erasing the node p'. Then we call (3 1 a reduction of f3 and call the way we obtained f3' 
from f3 a reduction step. A reduction procedure is a sequence of successive reduction 
steps. 

One sees now easily that a reduction procedure applied to (3 produces a new 
parameterized arithmetic circuit (3* (also called a reduction of /3) with the same basic 
parameter and input nodes, which computes the same final results as (3 (although 
their possible indexing may be changed). Moreover, if (3 is a robust parameterized 
circuit, then (3* is robust too. Observe also that in the case of robust parameterized 
circuits our reduction commutes with restriction. 

Now we introduce the second rewriting procedure. 
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Let assumptions and notations be as before and let be given a set P of nodes 
of and a robust parameterized arithmetic circuit 7 with parameter domain Ai 
and j^P input nodes, namely for each p G P one which becomes labelled by a new 
input variable Y p . We obtain a new parameterized arithmetic circuit, denoted by 
7 *p (3, when we join 7 with f3, replacing for each p G P the input node of 7, which 
is labelled by the variable Y p , by the node p of (3. The output nodes of (3 constitute 
also the output nodes of 7 *p (3. Thus (3 and 7 *p (3 compute the same final results. 
Observe that 7 *p [3 is robust if it is consistent. We call the circuit 7 *p [3 and all its 
reductions broadcastings of (3. Thus broadcasting a robust parameterized arithmetic 
circuit means rewriting it using only valid polynomial identities. 

If we consider arithmetic circuits as computer programs, then reduction and 
broadcasting represent a kind of program transformations. 



3.1.2 A specification language for circuits 

Computer programs (or "programmable algorithms" ) written in high level languages 
are not the same thing as just "algorithms" in Complexity Theory. Whereas in the 
uniform view algorithms become implemented by suitable machine models and in 
the non-uniform view by devices like circuits; specifications and correctness proofs 
are not treated by the general theory, but only, if necessary, outside of it in a case- 
by-case ad-hoc manner. The meaning of "algorithm" in Complexity Theory is 
therefore of syntactic nature. 

On the other hand, computer programs, as well as their subroutines (modules) 
include specifications and correctness proofs, typically written in languages orga- 
nized by a hierarchy of different abstraction levels. In this sense programmable 
algorithms become equipped with semantics. This is probably the main difference 
between Complexity Theory and Software Engineering. 

In this paper, we are only interested in algorithms which in some sense are pro- 
grammable. The routines of our computation model will operate on parameterized 



arithmetic circuits (see Section 3.3). Therefore we are now going to fix a (many- 



sorted) first-order specification language C for these circuits. 

The language C will include the following non-logical symbols: 

- 0, 1, +, — , x, and a constant for each complex number, 

- variables 
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x« ...,x^... 
r«...,yw... 

to denote non-negative integers and vectors of them, robust parameterized 
arithmetic circuits, their nodes, their parameter domains, their parameter in- 
stances, their input variable vectors and instances of input variable vectors in 
suitable affine spaces, 

- suitable binary predicate symbols to express relations like "p is a node of the 
circuit "multiplication is the label of the node p of the circuit U M 
is the parameter domain of the circuit (3" , U U is a parameter instance of the 
circuit /3" , "r is a non-negative integer and the vector length of U is r", U X 
is the input variable vector of the circuit f3" and "n is a non-negative integer 
and the vector length of X is n" , 

- a ternary predicate symbol to express "pi and P2 are nodes of the circuit (3 
and there is an edge of (3 from p\ to p 2 ", 

- binary function symbols to express U U is a parameter instance, A; is a natural 
number and Uk is the /c-th entry of U" and "X is an input variable vector, 
n is a natural number and X n is the n-th entry of X" and "Y is a variable 
vector instance, n is a natural number and Y n is the n-th entry of Y" , 

- a unary function and a binary predicate symbol to express "the set of output 
nodes of the circuit (3" and "p is an output node of the circuit (3" 

- a quaternary function symbol G p ((3; U; X) to express "p is a node of the circuit 
(3, U is a parameter instance and X is the input variable vector of (3 and 
G p ((3; U;X) is the intermediate result of (3 at the node p and the parameter 
instance [/", 

- a predicate symbol for equality for any of the sorts just introduced. 

For the treatment of non-negative integers we add the Presburger arithmetic to 
our first-order specification language C 

At our convenience we may add new function and predicate symbols and variable 
sorts to C Typical examples are for (3 a circuit, U a parameter instance, X the 
input variable vector and p, pi, . . . , p m nodes of (3: 

"degree of G p ((3\ U\X) n and "the vector lengths of X and Y are equal (say n) and 
Y is a point of the closed subvariety of A n defined by the polynomials G Pl (/3; U;X), 
...,G p J(3;U;Xy. 

In the same spirit, we may increase the expressive power of £ in order to be 
able to express for a robust parameterized circuit (3 with irreducible parameter 
domain, U a parameter instance, X the input variable vector, p a node of (3 and a 
a vector of non-negative integers of the same length as X (say n), "the coefficient 
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of the monomial X a occurring in the polynomial G P ((3;U; X)" (recall Lemma [7]). 
Here we denote for X := (Xi, . . . , X n ) and a := {ot\, . . . , a n ) by X a the monomial 

•V" ; x? v;;-. 

The semantics of the specification language £ is determined by the universe of all 
robust parameterized arithmetic circuits, where we interpret all variables, function 
symbols and predicates as explained before. We call this universe the standard 
model of C. The set of all closed formulas of C which are true in this model form 
the elementary theory of C. 



3.2 Generic computations 

In the sequel, we shall use ordinary arithmetic circuits over C as generic compu- 
tations [BCS97J (also called computation schemes in |Hei89j ). The indegree zero 
nodes of these arithmetic circuits are labelled by scalars and parameter and input 
variables. 

The aim is to represent different parameterized arithmetic circuits of similar size 
and appearance by different specializations (i.e., instantiations) of the parameter 
variables in one and the same generic computation. For a suitable specialization 
of the parameter variables, the original parameterized arithmetic circuit may then 
be recovered by an appropriate reduction process applied to the specialized generic 
computation. 

This alternative view of parameterized arithmetic circuits will be fundamental 
for the design of routines of the branching-free computation model we are going to 



describe in Section 3.3.2 The routines of our computation model will operate on 
robust parameterized arithmetic circuits and their basic ingredients will be subrou- 
tines which calculate parameter instances of suitable, by the model previously fixed, 
generic computations. These generic computations will be organized in finitely many 
families which will only depend on a constant number of discrete parameters. These 
discrete families constitute the basic building block of our model for branching-free 
computation. 

We shall now exemplify these abstract considerations in the concrete situation 
of the given parameterized arithmetic circuit /3. Mutatis mutandis we shall follow 
the exposition of |KP96j . Section 2. Let I, L , . . . , Li +1 with L > r + n + 1 and 
Li+\ > q be given natural numbers. Without loss of generality we may suppose that 
the non-scalar depth of f3 is positive and at most /, and that (3 has an oblivious 
levelled structure of I + 2 levels of width at most L Q , . . . , Li + i. Let Ui,...,U r be 
new indeterminates (they will play the role of a set of "special" parameter variables 
which will only be instantiated by tti, . . . , ir r ). 

We shall need the following indexed families of "scalar" parameter variables 
(which will only be instantiated by complex numbers): 

- for n + r < j < L the indeterminate V}; 
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- for 1 < i < I, 1 < j < Li, < h < i, 1 < k < L h , the indeterminates A\j k \ 
B«f and S id , V,,: 

- for 1 < j < Li +1 , 1 < k < Li the indeterminate Cj\ 

We consider now the following function Q which assigns to every pair (i,j), 
1 < % < I, 1 < j < Lj, and (/ + 1, j), 1 < j < L l+X the rational expressions defined 
below: 

Qo,i '■= Ui,..., Qo,r := U r , 

Qo,r+l '■— Xi, . . . , Qo,r+n '■— X n , 
Qo,r+n+l '■= K+n+1) • • • 5 Qo,L := ^L - 

For 1 < i < I and 1 < j < L; L the value Q it j of the function Q is recursively 
defined by 

Qi,j '■— Si,j( ^2 ^iJ^QKk ■ ^2 ^ij ' k ^Qh',k') + 

0<h<i 0<k'<i 

l<k<L h l<k'<L h , 

0<h<i 0<h'<i 

l<k<L h l<k'<L h , 

Finally, for (/ + 1, j), 1 < j < L l+1 we define Q(i + i d ) := Ei<fc<L ; CiQi.fc- 

We interpret the function Q as a (consistent) ordinary arithmetic circuit, say T, 
over Z (and hence over C) whose indegree zero nodes are labelled by the "standard" 
input variables Xi, . . . , X n , the special parameter variables U\, . . . , U r and the scalar 
parameter variables just introduced. 

We consider first the result of instantiating the scalar parameter variables con- 
tained in T by complex numbers. We call such an instantiation a specialization 
of T. It is determined by a point in a suitable afline space. Not all possible spe- 
cializations are consistent, giving rise to an assignment of a rational function of 
C(Ui, . . . ,U r ,Xi, . . . , X n ) to each node of T as intermediate result. 

We call the specializations which produce a failing assignment inconsistent. If 
in the context of a given specialization of the scalar parameter variables of T we 
instantiate for each index pair {i, j), 1 < % < I, 1 < j ' < Li the variables Sij and Tij 
by two different values from {0, 1}, the labelled directed acyclic graph T becomes 
an ordinary arithmetic circuit over C of non-scalar depth at most / and non-scalar 
size at most L ± + • • • + L t with the inputs U\, . . . , U r , Xi, . . . , X n . 

We may now find a suitable specialization of the circuit T into a new circuit T' 
over C such that the following condition is satisfied: 

the (by Ai) parameterized circuit obtained from V by replacing the special parame- 
ter variables Ui, . . . , U r by 7Ti, . . . , 7i>, is consistent and can be reduced to the circuit 

0- 
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We may consider the circuit T as a generic computation which allows to recover 
(3 by means of a suitable specialization of its scalar and special parameter vari- 
ables into complex numbers and basic parameters ttx, . . . , 7r r and by means of circuit 
reductions. Moreover, any by A4 parameterized, consistent arithmetic circuit of 
non-scalar depth at most I, with inputs X\, . . . ,X n and q outputs, which has an 
oblivious level structure with I + 2 levels of width at most Lq, . . . ,Li + \, may be 
recovered from T by suitable specializations and reductions (see |BCS97j . Chapter 
9 for more details on generic computations). 

3.3 A model for branching— free computation. 

3.3.1 Requirements to be satisfied by our branching— free computation 
model. Informal discussion. 

We are now going to introduce a model of branching-free computation with param- 
eterized arithmetic circuits. We shall first require that the routines of this compu- 
tation model should be well behaved under restrictions of the inputs. We discuss 
this issue first informally. 

Suppose for the moment that our branching-free computation model is already 
established. Then its routines transform a given robust parameterized arithmetic 
(input) circuit into another parameterized (output) circuit such that both circuits 
have the same parameter domain. Well behavedness under restrictions will be a 
property of circuit transformation that guarantees that the output circuit is still 
robust. In particular, we wish that the following requirement is satisfied. 

Let A be a routine of our branching-free computation model and consider the 
previously introduced parameterized circuit (3. Let Af be a constructible subset of 
AA and suppose that (3 is an admissible input for the routine A. Then A produces 
on input (3 a parameterized arithmetic output circuit with parameter domain AA 
which we denote by A{[3). In order to formulate for the routine A our requirement, 
we must be able to restrict (3 and A(/3) to Af. Thus (3 and A((3) should be robust, (3^ 
should be a consistent admissible input circuit for A and A({3jv) should be consistent 
too. 

Our architectural requirement on the routine A may now be formulated as fol- 
lows: 

The parameterized arithmetic circuit A{(3j^) can be recovered from A{(3) 
by restriction to Af and circuit reduction. 

Routines which are well behaved under restrictions will automatically satisfy this 
requirement. 

The routine A performs with the parameterized arithmetic circuit (3 a trans- 
formation whose crucial feature is that only nodes which depend on the inputs 
Xi, . . . ,X n of (3 become modified, whereas parameter nodes remain substantially 
preserved. This needs an explicitation. 
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Suppose that (3 has t essential parameter nodes. Then the essential parameters 
(intermediate results) of (3 associated with these nodes define a geometrically robust 
constructible map 9 : M. — > A*. The image T of 9 is a constructible subset of A*. 
We require now that, as far as A performs arithmetic operations with parameters 
of (3, A does it only with essential ones, and that all essential parameters of A(/3) 
are obtained in this way. Further we require that there exists a geometrically robust 
constructible map v defined on T (e.g., a polynomial map) such that the results of 
these arithmetic operations occur as entries of the composition map v o 9. From 
Corollary [6] we deduce that v o 9 is a geometrically robust constructible map. 

Our basic construction method of routines will be recursion. A routine of our 
computation model which can be obtained in this way is called recursive. 

Suppose now that A is a recursive routine of our computation model. Then A 
should be organized in such a way that for each internal node p of /3, which depends 
on at least one input, there exists a set of nodes of A((3), also denoted by p, with 
the following property: 

the elements of the set p of nodes of A((3) represent the outcome of the action of A 
at the node p of (3. 

We fix now a node p of (3 which depends on at least one input. Let G p be the 
intermediate result associated with the node p of (3 and let F p be a vector whose 
entries are the intermediate results of *4(/3) at the nodes contained in the set p 
of nodes of A{(3). Thus F p is a vector of rational functions in a suitable tuple of 
(standard) variables, say X' . 

Recall that by assumption (3 and A((3) are robust parameterized arithmetic cir- 
cuits with parameter domain M. . Therefore we deduce from Lemma [7] that G p and 
the entries of F p are in fact polynomials in X±, . . . ,X n and X', respectively, and 
that their coefficients are geometrically robust functions defined on A4. 

As part of our second and main requirement of our computation model we de- 
mand now that A satisfies at the node p of (3 the following isoparametricity condition: 

(i) for any two parameter instances U\ and w 2 of Ai the assumption 

G p (ui, X\, . . . , X n ) = Gp(w 2 , Xi, . . . , X n ) 

implies 

F p (u 1 ,X / ) = F p (u 2 ,X / ). 

Let 9 P be the coefficient vector of G p and observe that 9 P is a geometrically robust 
constructible map defined on Ml, whose image, say T p , is an irreducible constructible 
subset of a suitable affine space. 

Since the first-order theory of the algebraically closed field C admits quantifier 
elimination, one concludes easily that condition (i) is satisfied if and only if there 
exists a constructible map a p defined on T p such that the composition map a p o 9 p 
(which is also constructible) represents the coefficient vector of (all entries of) F p . 
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In the sequel we shall need that the dependence a p of the coefficient vector of 
F p on the coefficient vector of G p is in some stronger sense uniform (and not just 
constructible) . Therefore we include the following condition in our requirement: 

(ii) the constructible map a p is geometrically robust. 

The map o p is uniquely determined by condition (i). Moreover, the map a p 
depends on the (combinatorial) labelled DAG structure of (3 below the node p, 
but not directly on the basic parameters 7Ti,...,7r r . This is the essence of the 
isoparametric nature of conditions (i) and (ii). We shall therefore require that our 
recursive routine is isoparametric in this sense, i.e., that A satisfies conditions (i) 
and (ii) at any internal node p of (3 which depends at least on one input. 

Observe that the geometrically robust constructible map a p (which depends on 
j3 as well as on p) is not an artifact, but emerges naturally from the recursive con- 
struction of a circuit semantic within the paradigm of object-oriented programming. 
To explain this, let notations and assumptions be as before and suppose that A is a 
isoparametric recursive routine of our model and that we apply A to the robust pa- 
rameterized arithmetic circuit (3. Let p again be a node of f3 which depends at least 
on one input. Let u be a parameter instance of M. and denote by (3^ u \ G p u \ A(/3)^ 
and F p u ^ the instantiations of (3, G p , A((3) and F p at u (observe that the interme- 
diate results of (3^ and A(f3)^ are well defined although we do not require that 
these circuits are consistent). Then the intermediate results of A((3)^ contained in 
F p depend only on the intermediate result G p u ^ of (3^ and not on the parameter 
instance u itself. In this spirit we may consider the sets V p := {G p u ^ ; u G A4} and 
$ p := {Fp U ^ ; u G M.} as abstract data types and (3 and A((3) as syntactic descrip- 
tions of two abstraction functions which associate to any concrete object u G 
the abstract objects G p u>> and F p u \ respectively. The identity map idj^ : M. — > Ai 
induces now an abstract function |MeyOO| from T p to $ p , namely a p : Y p — > $ p . In 
this terminology, id M is just an implementation of a p . If we now consider that each 
recursive step of the routine A on input (3 has to be realized by some routine of the 
object-oriented programming paradigm, we arrive to a situation which requires the 
existence of a geometrically robust constructible map a p : Y p — >■ <3> p as above. 

We may interpret the map a p : Y p — >• $ p also as an ingredient of a specification of 
the recursive routine A. The map a p may be thought as an operational specification 
which determines F p in function of G p . A weaker specification would be a descriptive 
one which relates G p and F p without determining F p from G p completely. 

In order to motivate the requirement that the recursive routine A should be 
isoparametric, we shall consider the following condition for recursive routines which 
we call well behavedness under reductions. 



We only outline here this condition and leave the details until Section 3.3.2 
Suppose now that we apply a reduction procedure to the robust parameterized 
input circuit (3 producing thus another robust, by M. parameterized circuit (3* which 
computes the same final results as (3. Then the reduced circuit (3* should also be 
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an admissible input for the routine A. We call the recursive routine A well behaved 
under reductions if on input (3 it is possible to extend the given reduction procedure 
to the output circuit A((3) in such a way, that the extended reduction procedure, 
applied to A({3), reproduces the circuit A((3*). 

Obviously well behavedness under reductions limits the structure of A{f3). Later, 
in Section 3.3.2[ we shall see that, cum grano salis, any recursive routine, which is 



well behaved under restrictions and reductions, is necessarily isoparametric. Since 
well behavedness under restrictions and reductions are very natural quality at- 
tributes for routines which transform robust parameterized arithmetic circuits, the 
weaker requirement, namely that recursive routines should be isoparametric, turns 
out to be well motivated. 



In Section 3.3.2 we shall formally introduce our branching-free computation 
model. We postpone for then the precise definitions of the notions of well behaved- 
ness under restrictions and reductions. 

There exists a second reason to limit the recursive routines of our branching- 
free computation model to isoparametric ones. Isoparametric recursive routines 
have considerable advantages for program specification and verification by means of 



Hoare Logics (see |Apt81| ). We shall come back to this issue in Section 3.3.2.4 



3.3.2 The branching— free computation model 

The computation model we are going to introduce in this and the next subsection 
will be comprehensive enough to capture the essence of all known circuit based elim- 
ination algorithms in effective Algebraic Geometry and, mutatis mutandis, also of 
all other (linear algebra and truncated rewriting) elimination procedures (see Sec- 
tions 3.3.3, ij |Mor03j . |Mor05j . and the references cited therein, and for truncated 
rewriting methods especially [DFGS91J). The only algorithm from symbolic arith- 
metic circuit manipulation which will escape from our model is the Baur-Strassen 
gradient computation |BCS97j . Chapter 7.2. 

In the sequel we shall distinguish sharply between the notions of input variable 
and parameter and the corresponding categories of circuit nodes. 

Input variables, called "standard", will occur in parameterized arithmetic circuits 
and generic computations. The input variables of generic computations will appear 
subdivided in three sorts, namely as "parameter" , "argument" and "standard" input 
variables. 

The branching-free computation model we are going to introduce in this subsec- 
tion will assume different shapes, each shape being determined by a finite number 
of a priori given discrete (i.e., by tuples of natural numbers indexed) families of 
generic computations. The labels of the inputs of the ordinary arithmetic circuits 
which represent these generic computations will become subdivided into parameter, 
argument and standard input variables. We shall use the letters like U, U', U", . . . 
and W, W, W" to denote vectors of parameters, Y, Y', Y", . . . and Z, Z', Z" to de- 
note vectors of argument and X, X', X", ... to denote vectors of standard input 
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variables (see Section 3.2). 

We shall not write down explicitly the indexations of our generic computations 
by tuples of natural numbers. Generic computations will simply be distinguished 
by subscripts and superscripts, if necessary. 

Ordinary arithmetic circuits of the form 

R Xl (W i; XW), R X2 (W 2 ;X&), ... 
R> X1 {W V -XV)), R> X2 (W 2 ,;X^), ... 

represent a first type of a discrete family of generic computations (for each variable 
Xi,X 2 , . . . , X n we suppose to have at least one generic computation). Other types 
of families of generic computations are of the form 

R + (W;U,Y;X), R' + (W';U',Y';X'), R\{W";U" \Y";X") 

R.fiW',U,Y',X), R' (W', U', Y'\ X ), R"^W";U\Y";X") 

R add (W- Y } Z; X), R' add (W>; Y', Z>; X'), K dd (W"; Y", Z"; X") . . . 

R mult (W;Y,Z ] X), R' mult (W';Y',Z';X% K ult (W"; Y", Z"; X") ... 



and 



R div (W;Y,Z;X), R' div (W';Y' , Z'; X'), R% V (W";Y" , Z"; X") 



Here the subscripts refer to addition of, and multiplication or division by a parameter 
(or scalar) and to essential addition, multiplication and division. A final type of 
families of generic computations is of the form 

R(W;Y;X), R'(W; Y'\ X'), R"(W"; Y"\ X"), . . . 



We recall from Section 3.3.1 that the objects handled by the routines of any shape 
of our computation model will always be robust parameterized arithmetic circuits. 
The inputs of these circuits will only consist of standard variables. 

From now on we have in mind a previously fixed shape when we refer to the 
branching-free computation model we are going to introduce. We start with a given 
finite set of discrete families of generic computations which constitute a shape as 
described before. 



3.3.2.1 The notion of well behavedness under restrictions 

A fundamental issue is how we recursively transform a given input circuit into 
another one with the same parameter domain. During such a transformation we 
make an iterative use of previously fixed generic computations. On their turn these 
determine the corresponding recursive routine of our branching-free computation 
model. 

We consider again our input circuit j3. We suppose that we have already chosen 
for each node p, which depends at least on one of the input variables Xi, . . . , X n , a 
generic computation 

R%(W p ;X^), 
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and that this choice was made according to the label of p, namely JQ, 1 < i < n, 
or addition of, or multiplication or division by an essential parameter, or essential 
addition, multiplication or division. Here we suppose that U p is a single variable, 
whereas W p , Y p , Z p and 1^ may be arbitrary vectors of variables. 

Furthermore, we suppose that we have already precomputed for each node p 
of /3, which depends at least on one input, a vector w p of geometrically robust 
constructible functions defined on Ai. If p is an input node we assume that w p is a 
vector of complex numbers. Moreover, we assume that the length of w p equals the 
length of the variable vector W p . We call the entries of w p the parameters at the 
node p of the routine A applied to the input circuit (3. 

We are now going to develop the routine A step by step. The routine A takes 
over all computations of (3 which involve only parameter nodes, without modifying 
them. 

Consider an arbitrary internal node p of j3 which depends at least on one input. 
The node p has two ingoing edges which come from two other nodes of j3, say pi 
and pi- Suppose that the routine A, on input (3, has already computed two results, 
namely F pi and F p2 , corresponding to the nodes pi and p 2 . Suppose inductively that 
these results are vectors of polynomials depending on those standard input variables 
that occur in the vectors of the form X( p '\ where p' is any predecessor node of p. 
Furthermore, we assume that the coefficients of these polynomials constitute the 
entries of a geometrically robust, constructible map defined on M.. Finally we 
suppose that the lengths of the vectors F pi and Y p (or U p ) and F p2 and Z p coincide. 

The parameter vector w p of the routine A forms a geometrically robust, con- 
structible map defined on M, whose image we denote by JC p . Observe that K p is a 
constructible subset of the affine space of the same dimension as the length of the 
vectors w p and W p . Denote by n p the vector of the restrictions to /C p of the canon- 
ical projections of this affine space. We consider /C p as a new parameter domain 
with basic parameters n p . For the sake of simplicity we suppose that the node p is 
labelled by a multiplication. Thus the corresponding generic computation has the 
form 



R^{W p -U p ,Y p -X^) 



(1) 



or 



R ( :L(w p ,y p ,z p ,x^). 



(2) 



Let the specialized generic computations 



R^(k p ,U p ,Y p ,X^) and R ( ± t ( Kp ,Y p , Z p , X^) 
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be the by JC P parameterized arithmetic circuits obtained by substituting in the 
generic computations ([I]) and Q for the vector of parameter variables W p the basic 
parameters K p . At the node p we shall now make the following requirement on the 
routine A applied to the input circuit j3: 

(A) The by 1C P parameterized arithmetic circuit which corresponds to the current 
case, namely 

RW(k p ;U p ,Y p ;XM) 

or 

should be consistent and robust. 

Observe that the requirement (A) is automatically satisfied if all the generic 
computations of our shape are realized by totally division-free ordinary arithmetic 
circuits. 

Assume now that the routine A applied to the circuit (3 satisfies the requirement 
(A) at the node p of f3. 

Let us first suppose that the node p is labelled by a multiplication involving 
an essential parameter. Recall that in this case we assumed earlier that the length 
of the vector F Pl is one, that F Pl is an essential parameter of and that the vec- 
tors F p2 and Y p have the same length. Joining now with the generic computation 
B>y(Wp\ U p , Y p ; X^) at W p , U p and Y p the previous computations of w p , F pi and F p2 , 
we obtain a parameterized arithmetic circuit with parameter domain A4, whose final 
results are the entries of a vector which we denote by F p . 

Suppose now that the node p is labelled by an essential multiplication. Recall 
again that in this second case we assumed earlier the vectors F pi and Y p and F P2 
and Z p have the same length. Joining with the generic computation 

R { lu(W P ,Yp,Z p ;X^) 

at W p , Yp and Z p the previous computations of w p , F Pl and F P2 we obtain also a 
parameterized arithmetic circuit with parameter domain Ai, whose final results are 
the entries of a vector which we denote again by F p . 

One deduces easily from our assumptions on w p , F Pl and F P2 and from the re- 
quirement (A) in combination with Lemma [7] and Corollary |6j that in both cases 
the resulting parameterized arithmetic circuit is robust if it is consistent. The other 
possible labellings of the node p by arithmetic operations are treated similarly. In 
particular, in case that p is an input node labelled by the variable Xj, 1 < i < n, 
the requirement (A) implies that the ordinary arithmetic circuit R^_(w p ; X^) is 
consistent and robust and that all its intermediate results are polynomials in X^ 
over C (although R ( p.(w p ;X^) may contain divisions). 



In view of our comments in Section 3.3.1 . we call the recursive routine A (on in 



put f3) well behaved under restrictions if the requirement (A) is satisfied at any node 
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p of (3 which depends at least on one input and if joining the corresponding generic 
computation with w p , F Pl and F P2 produces a consistent circuit (observe that this 
last condition is automatically satisfied when the specialized generic computation of 
(A) is essentially division-free). If the routine A is well behaved under restrictions, 
then A transforms step by step the input circuit (3 into another consistent robust 
arithmetic circuit, namely *4(/3), with parameter domain Ai. Thus, well behaved- 
ness under restrictions guarantees that the recursive routine A transforms robust 
parameterized arithmetic circuits in robust ones. 

As a consequence of the recursive structure of A((3), each node p of (3 generates a 
subcircuit of A(/3) which we call the component of A((3) generated by p. The output 
nodes of each component of A((3) form the hypernodes of a hypergraph Ha(P) whose 
hyperedges are given by the paths connecting the nodes of A((3) contained in distinct 
hypernodes of W^iP)- The hypergraph Ha(/3) m8b Y be shrunk to the DAG structure 
of j3 and therefore we denote the hypernodes of 1~La(P) i n the same way as the nodes 
of f3. Notice that well behavedness under restrictions is in fact a property which 
concerns the hypergraph Ha(P)- 

We call A a (recursive) parameter routine if A does not introduce new standard 
variables. In the previous recursive construction of the routine A, the parameters 
at the nodes of (3, used for the realization of the circuit A((3), are supposed to be 
generated by recursive parameter routines. 



3.3.2.2 The notion of isoparametricity 

We are now going to consider another requirement of our recursive routine A, 
which will lead us to the notion of isoparametricity of A. Isoparametricity will 



guarantee that the recursive routine A may be specified (see Section 3.3.2.4). 

Let us turn back to the previous situation at the node p of the input circuit 
p. Notations and assumptions will be the same as before. From Lemma [7] we 
deduce that the intermediate result of (3 associated with the node p, say G p , is 
a polynomial in X\, . . . ,X n whose coefficients form the entries of a geometrically 
robust, constructible map defined on M., say 6 P . Let T p be the image of this map and 
observe that T p is a constructible subset of a suitable affine space. The intermediate 
results of the circuit A((3) at the elements of the hypernode p of %a(P) constitute a 
polynomial vector which we denote by F p . 

We shall now make another requirement on the routine A at the node p of (3. 

(B) There exists a geometrically robust constructible map a p defined on T p such 
that a p o 9 p constitutes the coefficient vector of F p . 

In view of the comments made in Section 13.3.11 we call the recursive routine A 
isoparametric (on input 0) if requirements (A) and (B) are satisfied at any node p 
of (3 which depends at least on one input. 

Let assumptions and notations be as before and consider again the node p of the 
circuit (3. Assume that the recursive routine A is well behaved under restrictions and 
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denote by r p the coefficient vector of F p . Observe that t p is a geometrically robust 
constructible map defined on M.. Assume, furthermore, that A, applied to the 
circuit 0, fulfils the requirement (B) at p. Then the topological robustness (which 
is a consequence of the geometrical robustness) of a p implies that the following 
condition is satisfied: 

{B') Let (uk)kem be a (not necessarily convergent) sequence of parameter instances 
Uk G M. and let u G M. such that {6 p {uk))km converges to 8 p {u). Then the 
sequence (r p (uk))keN is bounded. 

Suppose now that the recursive routine A is well behaved under restrictions and 
satisfies instead of (B) only condition (B') at the node p of (3. Let u G Ai be an 
arbitrary parameter instance. Then Theorem [3] implies that r p takes on the set 
{v! G M; 9 p {u') = p (u)} only finitely many values. In particular, for VJl u being the 
vanishing ideal of the C-algebra C[9 P ] at p (u), the entries of r p are integral over the 
local C-algebra C[0 p ]ot u (the argument for that relies on Zariski's Main Theorem 
and is exhibited in [C GH + 03] . Sections 3.2 and 5.1). This algebraic characterization 
implies that for given u G M. all the sequences {T p {uk))km of condition (B') have 
only finitely many distinct accumulation points. This shows that requirement (B) 
and condition (B') are closely related. 

Adopting the terminology of (GHMSllj, we call A coalescent (on input (3), if A 
is well behaved under restrictions and satisfies condition (B') for any node p of f3. 
Thus isoparametricity implies coalescence for A, but not vice versa. Nevertheless the 
notions of isoparametricity and coalescence become quite close for recursive routines 
which are well behaved under restrictions. 



3.3.2.3 The notion of well behavedness under reductions 

Suppose again that the recursive routine A is well behaved under restrictions. 
We call A well behaved under reductions (on input (3) if A((3) satisfies the following 
requirement: 

Let p and p' be distinct nodes of (3 which compute the same intermediate 
results. Then the intermediate results at the hypernodes p and p' ofH-Mp) 
are identical. Mutatis mutandis the same is true for the computation of 
the parameters of A at any node of (3. 

Assume that the routine A is recursive and well behaved under reductions. One 
verifies then easily that, taking into account the hypergraph structure %a{P) °f ^(13), 
any reduction procedure on (3 may canonically be extended to a reduction procedure 
of A(f3). 



In Section 3.3.1 we claimed that, cum grano salis, the requirement of well be- 
havedness under reductions implies the requirement of isoparametricity for recursive 
routines. We are going now to prove this. 
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Let notations and assumptions be as before and let us analyze what happens 
to the recursive routine A at the node p of f3. For this purpose we shall use the 
following broadcasting argument. 

Recall that G p and the entries of F p are the intermediate results of (3 and A((3) 
associated with p, where p is interpreted as a node of the input circuit (3 in the 
first case and as a hypernode of Ha(p) m the second one. Moreover recall that 
G p is a polynomial in Xi, . . . , X n , that the geometrically robust, constructible map 
9 p , defined on Ai, represents the coefficient vector of G p and that the irreducible 
constructible set T p is the image of 9 p . Observe that the entries of 9 p may be 
computed from 7Ti, . .. ,7T r by a robust arithmetic circuit (e.g., by interpolation of 
G p in sufficiently generic points of A n ). We consider now the robust parameterized 
arithmetic circuit 7 P which realizes the following trivial evaluation of the polynomial 
G p . 

- compute simultaneously from 7Ti, . . . , ir r all entries of 9 P and from X±, . . . , X n 
all monomials occurring in G p 

- compute Gp as a linear combination of the monomials of G p using as coefficients 
the entries of 9 p . 

The circuit 7 P has a single output node, say p', which computes the polynomial 

Gp. 

Now we paste, as disjointly as possible, the circuit 7 P to the circuit j3 obtaining 
thus a new robust, parameterized arithmetic circuit f3 p with parameter domain M.. 
Observe that (3 P contains (3 and j p as subcircuits and that p and p' are distinct nodes 
of f3 p which compute the same intermediate result, namely G p . The entries of 9 P are 
essential parameters of 7 P and hence also of f3 p . We suppose now that f3 p is, like (3, 
an admissible input for the recursive routine A. Let F p i be a vector whose entries 
are the intermediate results at the nodes of A(f3 p ) contained in the hypernode p' 
of Ha(p ). Analyzing now how A operates on the structure of the subcircuit 7 P of 
(3 P , we see immediately that there exists a geometrically robust constructible map 
o p defined on T p such that the composition map o p o 9 P constitutes the coefficient 
vector of F p i. Since by assumption the recursive routine A is well behaved under 
reductions and the intermediate results of (3 p at the nodes p and p' consist of the 
same polynomial G p , we conclude that the intermediate results at the hypernodes 
p and p' of 1~La{P ) are also the same. Therefore we may assume without loss of 
generality F p = F p /. Hence the geometrically robust, constructible map a p o 9 P 
constitutes the coefficient vector of F p . 

This proves that the recursive routine A satisfies, on input (3 and at the node 
p, the requirement (£>). Since j3 was an arbitrary admissible input circuit for the 
recursive routine A and p was an arbitrary node of (3 which depends on at least one 
input, we may conclude that A is isoparametric. The only assumption we made to 
draw this conclusion was that the extended circuit (3 P is an admissible input for the 
routine A. This conclusion is however not very restrictive because (3 and (3 p compute 
the same final results. 



29 



3.3.2.4 Isoparametricity and program specification 

In Section 3.3.1 , we mentioned that isoparametric routines are advantageous for 
program specification and verification. We are now going to explain this. 

Let notations and assumptions be as before and let in particular A be a recursive 
routine of our computation model which behaves well under restrictions. Assume 
that (3 is an admissible input for A and consider the specification language £ intro- 



duced in Section [3. 1.2[ Suppose that the routine A is given by an asserted program 
I! formulated in the elementary Hoare Logics of C (|Apt81|). The standard model 
of the elementary theory of C provides us with the states which define the semantics 
of n. The asserted program II represents the routine A as a loop which transforms 
node by node the labelled DAG structure of (3 into the labelled DAG structure of 

At each step of the loop a purely syntactic action, namely a graph manipulation, 
takes place. This action consists of the join of two or more labelled directed acyclic 
graphs. Simultaneously, in order to guarantee the correctness of the program IT, a 
loop invariant, formulated in our specification language £, has to be satisfied. 

This involves the semantics of C consisting of the universe of all robust pa- 
rameterized arithmetic circuits. A loop invariant as above is given by a formula 
f\(/3\, (3 2 , -Mi, Pi) of C containing the free variables /3i, (3 2 for circuits over the same 
parameter domain Ai\ and p\ for a node of (3i and a linked hypernode of (3 2 , such 
that these free variables become instantiated by (3, A((3), M. and the node p of (3 or 
the hypernode p of A(f3). The variables U^ l \ . . . , U^ m \ . . . and the standard input 
variable vectors . . . , X^ h \ . . . occur only bounded in /\([3i, (3%, Aii, pi) and the 
variables pi, . . . , pi, . . . occur all bounded except one, namely p%. 

For 7r := (tti, . . . , 7r r ) and given variables X, X' and p expressing a parameter 
instantiation, the input variable vectors of (3 and A{(3) and a node of f3, we denote 
by G p ((3] 7r; X) and F p (A((3); tt; X') the function symbols (or vectors of them) which 
express the intermediate results of (3 or A((3) corresponding to p. 

We require now that any formula of C built up by G pi , . . . , G Pl and , . . . , F p i i , 
and containing only (3, Ai and pi as free variables is equivalent to a formula built 
up only by G Pl , . . . , G Pl and G ^ , . . . , G p ^ . This implies that in C the intermediate 
result F p of A((3) is definable in terms of the intermediate result G p of (3. Applied 
to the node p of the concrete circuit j3 with parameter domain Ai, this means that 
for 6 P and r p being the coefficient vectors of G p ((3, it, X) and F p (A((3), ir, X') and T p 
being the image of 8 p , there exists a constructible map a p with domain of definition 
Tp such that r p = Op o 6 p holds. In particular, for u', u" G M. the assumption 
e p ( u ') = 6p{u") implies t p {u') = r p (u"). 

For the modelling of elimination algorithms this is a reasonable requirement (see 
Section [4|. If we require additionally that the transformation of G p ((3, ir, X) into 
Fp(A((3), ii, X') is continuous, then the constructible map a p has to be geometrically 
robust (see Section 3.3.1). 

In terms of the specification language C, this reasoning may be formulated as 
follows. 



30 



Let /3i,/?2, M-i and p\ be variables for robust parameterized arithmetic circuits, 
their parameter domains and their (hyper)nodes. We assume that there exist a 
formula 

in the free variables (3\, (32, Mi, pi such that for any concrete, for A admissible circuit 
(3 with parameter domain M. and basic parameter vector ir and for any node p of /3 
the following condition is satisfied: 

(*) £l(/3, A((3),M, p) determines the polynomial F P (A((3), ir, X') 
in terms of G p ((3, ir, X). 

If C and A satisfy this assumption we say in the spirit of Hoare Logics that C is 
expressive for the routine A. 

Observe that condition (*) guarantees that a postcondition for the circuit A({3) 
can always be translated into an equivalent precondition for the circuit 0. 

3.3.2.5 Operations with routines 

Let A and B be recursive routines as before and suppose that they are well 
behaved under restrictions and isoparametric or even well behaved under reductions. 
Assume that A{(3) is an admissible input for B. We define the composed routine 
B o A in such a way that (B o A){(3) becomes the parameterized arithmetic circuit 
B(A((3)). Since the routines A and B are well behaved under restrictions, we see 
easily that (B o A){(3) is a consistent, robust parameterized arithmetic circuit with 
parameter domain Ai . From Lemma [7] and Corollary [6] we deduce that B o A is a 
isoparametric recursive routine if A and B are isoparametric. In case that A and B 
are well behaved under reductions, one verifies immediately that B o A is also well 
behaved under reductions. Therefore, under these assumptions, we shall consider 
B o A also as a routine of our computation model. 

Unfortunately, the composition of two arbitrary coalescent recursive routines 
need not to be coalescent. Therefore we shall focus in the sequel our attention 
on isoparametric recursive routines as basic building blocks of the branching-free 
computation model we are going to introduce. 

The identity routine is trivially well behaved under restrictions and reductions 
and in particular isoparametric. 

Let A and B be two routines of our computation model and suppose for the sake 
of simplicity that they are recursive and well behaved under restrictions. Assume 
that the robust parameterized arithmetic circuit /3 is an admissible input for A and 
B and that there is given a one-to-one correspondence A which identifies the output 
nodes of A({3) with the input nodes of B((3). Often, for a given input circuit (3, the 
correspondence A is clear by the context. If we limit ourselves to input circuits j3 
where this occurs, we obtain from A and B a new routine, called their join, which 
transforms the input circuit (3 into the output circuit B((3) *\A((3) (here we suppose 
that B{(3) *\ A((3) is consistent). Analyzing now B((3) *\ A((3), we see that the join 
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of A with B is well behaved under restrictions in the most obvious sense. Since by 
assumption the routines A and B are recursive, the circuits A((3) and B((3) inherit 
from (3 a superstructure given by the hypergraphs T-La(P) an d "Hb(/3)- Analyzing 
again this situation, we see that any reduction procedure on (3 can be extended in 
a canonical way to the circuit B((3) *\ A(f3). This means that the join of A with 
B is also well behaved under reductions if the same is true for A and B. More 
caution is at order with the notions of isoparametricity and coalescence. In a strict 
sense, the join of two isoparametric or coalescent recursive routines A and B is not 
necessarily isoparametric or coalescent. However, conditions (B) or (B') are still 
satisfied between the output nodes of f3 and B((3) * A A(/3). A routine with one of 
these two properties is called output isoparametric or output coalescent, respectively. 

The union of the routines A and B assigns to the input circuit (3 the juxtaposition 
of A(f3) and B((3). Thus, on input (3, the final results of the union of A and B are the 
final results of A{j3) and B((3) (taken separately in case of ambiguity). The union 
of A and B behaves well under restrictions and reductions and is isoparametric if 
the same is true for A and B. 

Observe also that for a recursive routine A which behaves well under restrictions 
and reductions the following holds: let (3 be a robust parameterized arithmetic circuit 
that broadcasts to a circuit j3* and assume that (3 and j3* are admissible circuits for 
A. Then A((3) broadcasts to A(0*). 

3.3.2.6 End of the description of the branching free computation model 

From these considerations we conclude that routines, constructed as before by 
iterated applications of the operations isoparametric recursion, composition, join 
and union, are still, in a suitable sense, well behaved under restrictions and output 
isoparametric. If only recursive routines become involved that behave well under 
reductions, we may also allow broadcastings at the interface of two such operations. 

This remains true when we introduce, as we shall do now, in our computational 
model the following additional type of routine construction. 

Let /3 be the robust, parameterized circuit considered before, and let R(W; Y; X) 
be a generic computation belonging to our shape list. Let wp be a precomputed 
vector of geometrically robust constructible functions with domain of definition Ai 
and suppose that wp and W have the same vector length and that the entries of 
wp are the final results of an output isoparametric parameter routine applied to the 
circuit (3. Moreover suppose that the final results of (3 form a vector of the same 
length as Y. 

Let K, be the image of wp. Observe that /C is a constructible subset of the afline 
space which has the same dimension as the vector length of W. Denote by k the 
vector of the restrictions to /C of the canonical projections of this afline space. We 
denote by R(k; Y; X) the ordinary arithmetic circuit over C obtained by substituting 
in the generic computation R(W;Y;X) the vector of parameter variables W by k. 
We shall now make the following requirement: 
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(C) The ordinary arithmetic circuit R(k;Y; X) should be consistent and robust. 

Observe that requirement (C) is obsolete when R(W; Y; X) is a totally division- 
free ordinary arithmetic circuit. 

Suppose now that requirement (C) is satisfied. A new routine, say B, is obtained 
in the following way: on input /3 the routine B joins with the generic computation 
R(W; Y; X) at W and Y the previous computation of wp and the circuit (3. 

From Lemma[7]and Corollary [6] we deduce that the resulting parameterized arith- 
metic circuit B((3) has parameter domain M. and is robust if it is consistent. We 
shall therefore require that B(f3) is consistent (this condition is automatically satis- 
fied if R(k; Y\ X) is essentially division-free). One sees immediately that the routine 
B is well behaved under restrictions and reductions and is output isoparametric. 

From now on we shall always suppose that all our recursive routines are isopara- 
metric and that requirement (C) is satisfied when we apply this last type of routine 
construction. 

An elementary routine of our simplified branching-free computation model is 
finally obtained by the iterated application of all these construction patterns, in 
particular the last one, isoparametric recursion, composition, join and union. As 
far as only recursion becomes involved that is well behaved under reductions, we 
allow also broadcastings and reductions at the interface of two constructions. Of 
course, the identity routine belongs also to our model. The set of all these routines 
is therefore closed under these constructions and operations. 

We call an elementary routine essentially division-free if it admits as input only 
essentially division-free, robust parameterized arithmetic circuits and all specialized 
generic computations used to compose it are essentially division-free. The outputs 
of essentially division-free elementary routines are always essentially division-free 
robust circuits. The set of all essentially division-free elementary routines is also 
closed under the mentioned constructions and operations. 

We have seen that elementary routines are, in a suitable sense, well behaved 
under restrictions. In the following statement we formulate explicitly the property 
of an elementary routine to be output isoparametric. This will be fundamental in 
our subsequent complexity considerations. 

Proposition 8 Let A be an elementary routine of our branching-free computation 
model. Then A is output isoparametric. More explicitly, let (3 be a robust, parame- 
terized arithmetic circuit with parameter domain M.. Suppose that (3 is an admissible 
input for A. Let 8 be a geometrically robust, constructible map defined on M. such 
that 9 represents the coefficient vector of the final results of (3 and let T be the image 
of 9. Then T is a constructible subset of a suitable affine space and there exists a 
geometrically robust, constructible map a defined on T such that the composition 
map o~ o 9 represents the coefficient vector of the final results of A{(3). 

A complete proof of this proposition follows with some extra work from our 
previous argumentation and will be omitted here. In case that A is a recursive 
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routine, Proposition [8] expresses nothing but the requirement (B) applied to the 
output nodes of 0. 

Let assumptions and notations be as in Proposition [8] and suppose that there is 
given a (not necessarily convergent) sequence (u^ken °f parameter instances u k G A4 
and that there exists a (possibly unknown) parameter instance u G Ai such that 
the sequence {6{uk))km converges to 6{u). In the spirit of |Ald84] . |Lic90] . §A and 
[BCS97J the sequence of (not necessarily consistent) ordinary arithmetic circuits 
(/^"*^)fceN represents an approximative algorithm for the instantiation of the final 
results of (3 at u. From Theorem [4] we conclude that the constructible map o 
is strongly continuous and therefore the sequence (A((3Y UkS) ) keN represents also an 
approximative algorithm for the instantiation of the final results of A((3) at u. 

One sees easily that this property characterizes output parametricity of routines 
which are well behaved under restrictions. 

Let us observe that Proposition [8] implies the following result. 

Corollary 9 Let assumptions and notations be as in Proposition^ Then the rou- 
tine A is output coalescent and satisfies the following condition: 

(*) Let u be an arbitrary parameter instance of M. and let VJt u be the vanishing 
ideal of the C-algebra C[9] at the point 9(u). Then the entries of the coefficient 
vector of the final results of A((3) are integral over the local C-algebra C[6]<xn u . 

The output coalescence of A and condition (*) are straight-forward consequences 
of the output isoparametricity of A. We remark here that condition (*) follows al- 
ready directly from the output coalescence of A. This highlights again the close con- 
nection between isoparametricity and coalescence. The argument requires Zariski's 
Main Theorem. For details we refer to [CG H + 03] . Sections 3.2 and 5.1. 

3.3.3 The extended computation model 

We are now going to extend our simplified branching-free computation model of 
elementary routines by a new model consisting of algorithms and procedures which 
may contain some limited branchings. Our description of this model will be rather 
informal. An algorithm will be a dynamic DAG of elementary routines which will be 
interpreted as pipes. At the end points of the pipes, decisions may be taken which 
depend on testing the validity of suitable universally quantified Boolean combina- 
tions of equalities between robust constructible functions defined on the parameter 
domain under consideration. The output of such an equality test is a bit vector 
which determines the next elementary routine (i.e., pipe) to be applied to the out- 
put circuit produced by the preceding elementary routine (pipe). This gives rise to a 
extended computation model which contains branchings. These branchings depend 
on a limited type of decisions at the level of the underlying abstract data type, 
namely the mentioned equality tests. We need to include this type of branchings in 
our extended computation model in order to capture the whole spectrum of known 
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elimination procedures in effective Algebraic Geometry. Because of this limitation 
of branchings, we shall call the algorithms of our model branching parsimonious 
(compare jGHOl] and |CGH+03j ). A branching parsimonious algorithm A which 
accepts a robust parameterized arithmetic circuit with parameter domain hA as 
input produces a new robust circuit A(0) with parameter domain M.. In particular 
A(/3) does not contain any branchings. 

Recall that our two main constructions of elementary routines depend on a pre- 
vious selection of generic computations from a given shape list. This selection may 
be handled by calculations with the indexing of the shape list. We shall think that 
these calculations become realized by deterministic Turing machines. At the begin- 
ning, for a given robust parametric input circuit with parameter domain Ai, a 
tuple of fixed (i.e., of independent) length of natural numbers is determined. This 
tuple constitutes an initial configuration of a Turing machine computation which de- 
termines the generic computations of our shape list that intervene in the elementary 
routine under construction. The entries of this tuple of natural numbers are called 
invariants of the circuit 0. These invariants, whose values may also be Boolean (i.e., 
realized by the natural numbers or 1), depend mainly on algebraic or geometric 
properties of the final results of 0. However, they may also depend on structural 
properties of the labelled DAG 0. 

For example, the invariants of may express that has r parameters, n inputs 
and outputs, (over C) non-scalar size and depth at most L and /, that is totally 
division-free, that the final results of have degree at most d < 2 l and that for 
any parameter instance their specializations form a reduced regular sequence in 
C[Xl, . . . , X n ], where Xi, . . . , X n are the inputs of 0. 

Some of these invariants (e.g., the syntactical ones like number of parameters, 
inputs and outputs and non-scalar size and depth) may simply be read-off from the 
labelled DAG structure of 0. Others, like the truth value of the statement that the 
specializations of the final results of at any parameter instance form a reduced 
regular sequence, have to be precomputed by an elimination algorithm from a previ- 
ously given software library in effective commutative algebra or algebraic geometry 
or their value has to be fixed in advance as a precondition for the elementary routine 
which becomes applied to 0. 

In the same vein we may equip any elementary routine A with a Turing com- 
putable function which from the values of the invariants of a given input circuit 
decides whether is admissible for A, and, if this is the case, determines the generic 
computations of our shape list which intervene in the application of A to 0. 

We shall now go a step further letting depend the internal structure of the com- 
putation on the circuit 0. In the simplest case this means that we admit that the 
vector of invariants of 0, denoted by inv(/3), determines the architecture of a first 
elementary routine, say A^m), which admits as input. Observe that the archi- 
tectures of the elementary routines of our computation model may be characterized 
by tuples of fixed length of natural numbers. We consider this characterization as 
an indexing of the elementary routines of our computation model. We may now use 
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this indexing in order to combine dynamically elementary routines by composition, 
join and union. Let us restrict our attention to the case of composition. In this 
case the output circuit of one elementary routine is the input for the next routine. 
The elementary routines which compose this display become implemented as pipes 
which start with a robust input circuit and end with a robust output circuit. Given 
such a pipe and an input circuit 7 for the elementary routine B representing the 
pipe, we may apply suitable equality tests to the final results of B(-y) in order to 
determine a bit vector which we use to compute the index of the next elementary 
routine (seen as a new pipe) which will be applied to -6(7) as input. 

A low level program of our extended computation model is now a text, namely 
the transition table of a deterministic Turing machine, which computes a function 
ip realizing the following tasks. 

Let as before (3 be a robust parameterized arithmetic circuit. Then if> returns 
first on input inv(/3) a Boolean value, zero or one, where one is interpreted as the 
informal statement u (3 is an admissible input". If this is the case, then if> returns 
on inv(/3) the index of an elementary routine, say A mv (/3), which admits (3 as input. 
Then ip determines the equality tests which have to be realized with the final results 
of Ai nv (p)(/3). Depending on the outcome of these equality tests ip determines an 
index value corresponding to a new elementary routine which admits -4 in v(/3)(/3) as 
input. Continuing in this way one obtains as end result an elementary routine 
A^\ which applied to f3, produces a final output circuit A^'(/3). The function ip 
represents all these index computations. We denote by if)(/3) the dynamic vector of 
all data computed by ip on input 0. 

The algorithm represented by ip is the partial map between robust parametric 
arithmetic circuits that assigns to each admissible input (3 the circuit A^\(3) as 
output. Observe that elementary routines are particular algorithms. This kind 
of algorithms constitute our extended computation model. We remark that any 
algorithm of this model is output isoparametric. If the pipes of an algorithm are all 
represented by essentially division-free elementary routines, we call the algorithm 
itself essentially division-free. 

One sees easily that the "Kronecker algorithm" |GLS01j (compare also |GHM + 98] , 
[GH H+97| and |GHMP97| 1 for solving non-degenerate polynomial equation systems 
over the complex numbers may be programmed in our extended computation model. 
Observe that the Kronecker algorithm requires more than a single elementary rou- 
tine for its design. In order to understand this, recall that the Kronecker algorithm 
accepts as input an ordinary division-free arithmetic circuit which represents by 
its output nodes a reduced regular sequence of polynomials G\, . . . ,G n belonging to 
C[Xi, . . . , X n ]. In their turn, the polynomials G\, . . . ,G n determine a degree pattern, 
say A := (Si, . . . , S n ), with Si := deg{Gi = 0, . . . , Gi — 0} for 1 < i < n. 

After putting the variables Xi, . . . , X n in generic position with respect to G\, . . . , 
G n , the algorithm performs n recursive steps to eliminate them, one after the other. 
Finally the Kronecker algorithm produces an ordinary arithmetic circuit which com- 
putes the coefficients of n + 1 univariate polynomials P,V\, . . . ,V n over C. These 
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polynomials constitute a "geometric solution" (see |GLS01j ) of the equation system 
G\ = 0, . . . , G n = because they represent the zero dimensional algebraic variety 
V := {G\ = 0, . . . , G n = 0} in the following "parameterized" form: 

V:={(V 1 (t),...,V n (t));teC,P(t) = 0}. 

Let (3 be any robust, parameterized arithmetic circuit with the same number of 
inputs and outputs, say X%, . . . , X n and Gx(U,Xi, . . . , X n ), . . . , G n (U, X%, . . . , X n ), 
respectively. Suppose that the parameter domain of (3, say M., is irreducible and 
that inv(/3) expresses that for each parameter instance u G Ai the polynomi- 
als Gi(u, Xi, . . . , X n ), . . . , G n (u, Xi, . . . , X n ) form a reduced regular sequence in 
C[Xi, . . . ,X n ] with fixed (i.e., from u G M. independent) degree pattern. Suppose, 
furthermore, that the degrees of the individual polynomials G\(u, Xi, . . . , X n ), . . . , 
G n (u,Xi, . . . ,X n ) are also fixed and that the variables Xi, . . . ,X n are in generic 
position with respect to the varieties {G\(u, X) = 0, . . . ,Gi(u,X) = 0}, 1 < i < n. 
Then, on input (3, the Kronecker algorithm runs a certain number (which depends on 
A) of elementary routines of our computation model which finally become combined 
by consistent iterative joins until the desired output is produced. 

Another non-trivial example for an algorithm of our extended computation 
model, which involves only limited branchings, is the Gaussian elimination pro- 
cedure of |Edm67j (or [Bar68j) applied to matrices whose entries are polynomials 
represented by ordinary arithmetic circuits in combination with a identity-to-zero 
test for such polynomials. The variables of these polynomials are considered as ba- 
sic parameters and any admissible input circuit has to satisfy a certain precondition 
formulated as the non-vanishing of suitable minors of the given polynomial matrix. 
Details and applications of this type of Gaussian elimination for polynomial matrices 
can be found in [Hei83j . 

3.3.3.1 Procedures 

We say that a given algorithm A of our extended model computes (only) param- 
eters if A satisfies the following condition: 

for any admissible input (3 the final results of A(j3) are all parameters. 

Suppose that A is such an algorithm and (3 is the robust parametric arithmetic 
circuit with parameter domain M. which we have considered before. Observe that 
A{(3) contains the input variables X\, . . . , X n and that possibly new variables, which 
we call auxiliary, become introduced during the execution of the algorithm A on in- 
put (3. Since the algorithm A computes only parameters, the input and auxiliary 
variables become finally eliminated by the application of recursive parameter rou- 
tines and evaluations. We may therefore collect garbage in order to reduce A(f3) to 
a final output circuit *4fmai(/3) whose intermediate results are only parameters. 

If we consider the algorithm A as a partial map which assigns to each admissible 
input circuit (3 its final output circuit -4fi na i(/3), we call A a procedure. 
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In this case, if ip is a low level program defining A, we call if) a low level procedure 
program. 

A particular feature of our extended computation model is the following: 
there exists a non-negative integer / (depending on the recursion depth of A) and 
non-decreasing real valued functions C/>0,...,Co>0 depending on one and the 
same dynamic integer vector, such that with the previous notations and Lp, L^p) 
denoting the non-scalar sizes of the circuits (3 and A{0) the condition 

L AW) <C f (4>({3))L f p + --- + C ^({3)) 

is satisfied. 

In the case of the Kronecker algorithm (and most other elimination algorithms 
of effective Algebraic Geometry) we have / := 1, because the recursion depth of the 
basic routines which intervene is one. 

In the sequel we shall need a particular variant of the notion of a procedure 
which enables us to capture the following situation. 

Suppose we have to find a computational solution for a formally specified gen- 
eral algorithmic problem and that the formulation of the problem depends on cer- 
tain parameter variables, say Ui, . . . , U r , input variables, say Xi, . . . , X n and output 
variables, say Yi, . . . , Y s . Let such a problem formulation be given and suppose that 
its input is implemented by the robust parameterized arithmetic circuit consid- 
ered before, interpreting the parameter variables Ui, . . . , U r as the basic parameters 

7Ti, . . . , 7T r . 

Then an algorithm A of our extended computation model which solves the given 
algorithmic problem should satisfy the architectural requirement we are going to 
describe now. 

The algorithm A should be the composition of two subalgorithms A^ 1 ' and A^ 
of our computation model which satisfy on input (3 the following conditions: 

(i) The subalgorithm A^ computes only parameters, (3 is admissible for A^ and 
none of the indeterminates Y\,...,Y S is introduced in A^\/3) as auxiliary 
variable (all other auxiliary variables become eliminated during the execution 
of the subalgorithm A^ on the input circuit (3). 

(ii) The circuit A^ al (f3) is an admissible input for the subalgorithm A^ 2 \ the in- 
determinates Yi, . . . , Y s occur as auxiliary variables in A^ {A^l al {f3)) and the 
final results of A^(Afl^ a i(0)) depend only on ttx, . . . , ir r and Y i: . . . , Y s . 

To the circuit A^^A^^P)) we may, as in the case when we compute only pa- 
rameters, apply garbage collection. In this manner A^^A^^)) becomes reduced 
to a final output circuit Aa na \(/3) with parameter domain M. which contains only 
the inputs Y±, . . . ,Y S . 

Observe that the subalgorithm A^ is by Proposition [8] an output isoparamet- 
ric procedure of our extended computation model (the same is also true for the 
subalgorithm A^ 2 \ but this will not be relevant in the sequel). 
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We consider the algorithm A, as well as the subalgorithms A^ and A^ 2 \ as 
procedures of our extended computation model. In case that the subprocedures 
A^ and A^ are essentially division-free, we call also the procedure A essentially 
division-free. This will be of importance in Section |4} 

The architectural requirement given by conditions (i) and (ii) may be interpreted 
as follows: 

the subprocedure A^ is a pipeline which transmits only parameters to the sub- 
procedure A^ . In particular, no (true) polynomial is transmitted from A^ to 

A®. 

Nevertheless, let us observe that on input (5 the procedure A establishes by 
means of the underlying low level program ip an additional link between j3 and the 
subprocedure A^ applied to the input A^(f3). The elementary routines which 
constitute A^ on input A^^fl) become determined by index computations which 
realizes ip on inv(/3) and certain equality tests between the intermediate results of 
A^\f3). In this sense the subprocedure A^ transmits not only parameters to the 
subprocedure but also a limited amount of digital information which stems from the 
input circuit (3. 

The decomposition of the procedure A into two subprocedures A^ and A^ sat- 
isfying conditions (i) and (ii) represents an architectural restriction which is justified 
when it makes sense to require that on input (3 the number of essential additions 
and multiplications contained in *4fi na i(/3) is bounded by a function which depends 
only on inv(/3). In Section 4.1, we shall make a substantial use of this restriction 
and give such a justification in the particular case of elimination algorithms. 

Here, we shall only point out the following consequence of this restriction. Let 
assumptions and notations be as before, let G, v and F be vectors composed by 
the final results of /3, A^\/3) and ^4fi na i(/3), respectively, and let 9 and (p be the 
coefficient vectors of G and F. Then the images of 9 and v are constructible subsets 
T and T 7 of suitable affine spaces and there exist geometrically robust constructible 
maps a and a 1 defined on T and T' with v = a o 9 and (p = a' o v = a' o a o 9. 

Based on [HK04] and |GHKaj . we shall develop in future work a high level 
specification language for algorithms and procedures of our computation model. 
The idea is to use a generalized variant of the extended constraint data base model 
introduced in [HK04] in order to specify algorithmic problems in symbolic Scientific 
Computing, especially in effective Algebraic Geometry. In this sense the procedure 
A, which solves the algorithmic problem considered before, will turn out to be 
a query computation composed by two subprocedures namely A^ and A^ which 
compute each a subquery of the query which specifies the given algorithmic problem. 
All these queries are called geometric because the procedures A^\ A^ and A are 
output isoparametric (see [ GHKaj ). 



39 



4 Applications of the extended computation model to 
complexity issues of effective elimination theory 

In this section we shall always work with procedures of our extended, branching 
parsimonious computation model. We shall study representative examples of elim- 
ination problems in effective Algebraic Geometry which certify, to a different ex- 
tent, that branching parsimonious elimination procedures based on our computation 
paradigm cannot run in polynomial time. 

4.1 Flat families of zero— dimensional elimination problems 

We first introduce, in terms of abstract data types, the notion of a flat family of 
zero-dimensional elimination problems (see also |GH01j and [CGH + 03J). Then we 
fix the classes of (concrete) objects, namely robust parameterized arithmetic circuits 
with suitable parameter domains, which represent ( "implement" ) these problems by 
means of a suitable abstraction function. 

Throughout this section, we suppose that there are given indeterminates U\, . . . , 
U r , X\, . . . , X n and Y over C. 

As concrete objects we shall consider robust parameterized arithmetic input and 
output circuits with parameter domain A r . The indeterminates Ui, . . . , U r will play 
the role of the basic parameters. The input nodes of the input circuits will be 
labelled by Xi, . . . ,X n , whereas the output circuits will have a single input node, 
labelled by Y. 

Let us now define the meaning of the term "flat family of zero-dimensional elim- 
ination problems" (in the basic parameters U\, . . . , U r and the inputs X 1; . . . , X n ). 
Let U := (Ui, . . . , U r ) and X := (X%, . . . , X n ) and let Gi, . . . , G n and H be poly- 
nomials belonging to the C-algebra C[U, X] := C[Z7i, . . . , U r , Xi, . . . , X n ]. Suppose 
that the polynomials G\, . . . ,G n form a regular sequence in C[U, X], thus defining an 
equidimensional subvariety V := {G\ = 0, . . . , G„ = 0} of the (n + r)-dimensional 
affine space A r x A n . The algebraic variety V has dimension r. Let 5 be the (geomet- 
ric) degree of V (observe that this degree does not take into account multiplicities or 
components at infinity). Suppose, furthermore, that the morphism of affine varieties 
7i : V — > A' r , induced by the canonical projection of A r x A n onto A r , is finite and 
generically unramified (this implies that ti is flat and that the ideal generated by 
G\, . . . , G n in C[U, X] is radical). Let t\ : V — > A r+1 be the morphism defined by 
tt(v) := (tt(u), H{y)) for any point v of the variety V. The image of tt is a hypersur- 
face of A r+1 whose minimal equation is a polynomial of C[U, Y] := C[Z7i, . . . , U r , Y] 
which we denote by F. Let us write deg F for the total degree of the polynomial F 
and degy F for its partial degree in the variable Y . Observe that F is monic in Y 
and that degF < 5 deg H holds. Furthermore, for a Zariski dense set of points u of 
A r , we have that deg y F is the cardinality of the image of the restriction of H to 
the finite set tt~ 1 (u). The polynomial F(U,H) vanishes on the variety V. 

Let us consider an arbitrary point u := (ui, . . . ,u r ) of A r . For given poly- 
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nomials A E C[U,X] and B E C[U,Y], we denote by and the poly- 
nomials A(u±, . . . , u r , Xi, . . . , X n ) and B(m, . . . , u r , Y) which belong to C[X] := 
C[Xi, . . . ,X n ] and C[Y] respectively. Similarly we denote for an arbitrary polyno- 
mial C E C[U] by the value C(ui, . . . ,u r ) which belongs to the field C. The 
polynomials G^\ . . . , G^ define the zero-dimensional subvariety 

V (u) ._ | G H = o, . . . , Gi M) = 0} ^ tT 1 ^) 

of the affine space A n . The degree (i.e., the cardinality) of V^ u ' is bounded by S. 
Denote by tp") : V^ u ' — > A 1 the morphism induced by the polynomial H^ u ' on the 
variety V^. Observe that the polynomial F^ u ' vanishes on the (finite) image of the 
morphism tt^. Observe also that the polynomial is not necessarily the minimal 
equation of the image of n^ u \ 

We call the equation system Gx — 0, . . . , G n — and the polynomial H a flat 
family of zero-dimensional elimination problems depending on the basic parameters 
Ui, . . . , U r and the inputs Xi, . . . , X n and we call F the associated elimination poly- 
nomial. A point u E A r is considered as a parameter instance which determines a 
particular problem instance, consisting of the equations G± = 0, . . . , G^ = and 
the polynomial H^ w \ A power of the polynomial F^ u ' is called a solution of this 
particular problem instance. 

We suppose now that the given flat family of elimination problems is imple- 
mented by an essentially division-free, robust parameterized arithmetic circuit (5 
with parameter domain A r and inputs X\, . . . , X n , whose final results are the poly- 
nomials Gi, . . . ,G n and H. The task is to find another essentially division-free, 
robust parameterized arithmetic circuit 7 with parameter domain A r having a sin- 
gle output node, labelled by Y, which computes for a suitable integer q E N the 
power F q of the associated elimination polynomial F. We suppose, furthermore, 
that this goal is achieved by the application of an essentially division-free procedure 
A of our extended computation model to the input circuit (3. Thus we may put 
7 := ^4finai(/3) and 7 may be interpreted as an essentially division-free circuit over 
C[U] with a single input Y (observe that the parameters computed by the robust 
circuits (3, A((3) and A$nai(P) are geometrically robust constructible functions with 
domain of definition A r which belong by [GHMSllJ, Corollary 12 to the C-algebra 
C [£/]). Using the geometric properties of flat families of zero-dimensional problems, 
we deduce from |GHM+98j . |GHH+97j . [GHMP97j . [GLSOlj or alternatively from 
|CGH89j . |DFGS91J that such essentially division-free procedures always exist and 
that they compute even the elimination polynomial F (the reader may notice that 
one needs for this argument the full expressivity of our computation model which 
includes divisions by parameters). 

We say that the essentially division-free procedure A solves algorithmically the 
general instance of the given flat family of zero-dimensional elimination problems if 
A computes F or a power of it. 

From now on we suppose that there is given a procedure A of our extended 
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computation model, d ecomp osed in two essentially division-free subprocedures A^ 
and A^ as in Section 3.3.3[ such that A solves algorithmically the general instance 
of any given flat family of zero-dimensional elimination problems. Our circuit (3 is 
therefore an admissible input for A and hence for A^. The final results of A^ l \f3) 
constitute a geometrically robust constructible map v defined on A r which repre- 
sents by means of A^ aY ((3) an admissible input for the procedure A^. Moreover, 
7 := A& na \((3) is an essentially division-free parameterized arithmetic circuit with 
parameter domain A r and input Y. 

Let S be the image of the geometrically robust constructible map v. Then S is an 
irreducible constructible subset of a suitable affine space. Analyzing now the internal 
structure of the essentially division-free, robust parameterized arithmetic circuit 
A^ (A^\(3)) , one sees easily that there exists a geometrically robust constructible 
map if) defined on S such that the entries of the geometrically robust composition 
map v* := if) o v constitute the essential parameters of the circuit 7. Let m be the 
number of components of the map v* . Since v and v* are composed by geometrically 
robust constructible functions defined on A r , we deduce from (11 IMS 1 1 j. Corollary 
12 that v and u* may be interpreted as vectors of polynomials of C[U]. 

The circuit 7 is essentially division-free. Hence there exists a vector u of m- 
variate polynomials over C such that the polynomials of C[Z7], which constitute 
the entries of u)(v*), become the coefficients of the elimination polynomial F with 
respect to the main indeterminate Y (see |KP96j . Section 2.1). Observe that we 
may write uj(u*) = u o v* interpreting the entries of v* as polynomials of C[U}. 

We are now going to see what happens at a particular parameter instance u G 
A r . Since (3, A^\(3), A(f3) and 7 = ^4fi na i(/3) are essentially division-free, robust 
parameterized arithmetic circuits with parameter domain A r , we may specialize the 
vector U of basic parameters to the parameter instance u G A r , obtaining thus 
ordinary division-free arithmetic circuits over C with the same inputs. We denote 
them by the superscript u, namely by (5 {u \ (A^ (j3))( u \ (^4(/3))( u ^ and 7^. One 
sees immediately that G[ u) , . . . , G { n u) and #(«> are the final results of (3^ u \ that the 
entries of v(u) are the final results of (A^ and that (F^) 9 is the final result 
of A{f3Y u " > and 7^. Observe that the division-free circuit 7 < - n - ) uses only the entries 
of u*(u) and fixed rational numbers as scalars. 

In the same spirit as before, we say that the procedure A solves algorithmically 
the particular instance, which is determined by u, of the given flat family of zero- 
dimensional elimination problems. 

Let us here clarify how all this is linked to the rest of the terminology used 
in |CGH + 03] . In this terminology the polynomial map given by to defines a "holo- 
morphic encoding" of the set of solutions of all particular problem instances and 
v*{u) is a "code" of the particular solution (F^) q . In the same context the ro- 
bust constructible map u* is called an "elimination procedure" which is "robust" 
since the procedure A^ is output isoparametric and since v* is geometrically robust 
(compare |CGH + 03"] . Definition 5, taking into account Lemma [7J Proposition [8] and 
Corollary [9] above) . 
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In this sense, we speak about families of zero-dimensional elimination problems 
and their instances and not simply about a single (particular or general) zero- 
dimensional elimination problem. 

Let us now turn back to the discussion of the given essentially division-free 
procedure A which solves algorithmically the general instance of any flat family of 
zero-dimensional elimination problems. 

We are now going to show the main result of this paper, namely that the given 
procedure A cannot run in polynomial time. 

Theorem 10 Let notations and assumptions be as before. For any natural num- 
ber n there exists an essentially division- free, robust parameterized arithmetic cir- 
cuit j3 n with basic parameters T, Ui, . . . , U n and inputs Xi, . . . ,X n which for U := 
(U\, . . . , U n ) and X := (X\, . . . ,X n ) computes polynomials G\ , . . . , Gn G C[X] 
and H^ n > G C[T, U,X] such that the following conditions are satisfied: 

(i) The equation system G± = 0, . . . , Gn = and the polynomial constitute 
a flat family of zero-dimensional elimination problems, depending on the pa- 
rameters T, Ui, . . . ,U n and the inputs X±, . . . , X n , with associated elimination 
polynomial G C[T, U, Y}. 

(ii) P n is an ordinary division-free arithmetic circuit of size 0(n) over C with 
inputs T, U 1 ,..., U n , Xi, ...,X n . 

(Hi) 7„ := Afi na i(f3 n ) is an essentially division-free, robust parameterized arithmetic 
circuit with basic parameters T, Ui, . . . , U n and input Y such that 7 n computes 
for a suitable integer g n eN the polynomial (F ( - n ^) qn . The circuit 7„ performs 
at least f2(2a ) essential multiplications and at least Q(2 n ) multiplications with 
parameters. Therefore ^ n has, as ordinary arithmetic circuit over C with in- 
puts T, Ui, . . . , U n and Y , non-scalar size at least Q(2 n ). 

Proof. During our argumentation we shall tacitly adapt to the new context the 
notations introduced before. We shall follow the main technical ideas behind the 
papers |GH01] . [CGH+03] and jGHMSll] . We fix now the natural number n and 
consider the polynomials 

c 1 ■ /^( n ) v n ■ r^i 71 ) ■ y 

Lri .— Lr 1 .— A 1 — Ai, . . . , Or n .— Or,^ .— Aj — A n 

and 

H := := T ~ 1X * + T II t 1 + & ~ X ) X i) 

l<i<n l<j<n 

which belong to C[X] and to C[T, U,X], respectively. 

Observe that Gi, . . . , G n and H may be evaluated by a division-free ordinary 
arithmetic circuit /3 := f3 n over C which has non-scalar size 0(n) and inputs T, 
U%, . . . , U n , Xi, . . . ,X n . As parameterized arithmetic circuit (3 is therefore robust. 
Hence satisfies condition (ii) of the theorem. 



43 



One sees easily that G\ — 0, . . . , G n — and H constitute a flat family of zero- 
dimensional elimination problems depending on the parameters T, U±, . . . ,U n and 
the inputs Xi, . . . , X n . 

Let us write H as a polynomial in the main indeterminates Xi,...,X n with 
coefficients Kli ..., Kn £ C[T, [/], «i, . . . , « n G {0, 1}, namely 

E ^,..,^r,---,^"- 

Ki,...,K n e{o,i} 

Observe that for k\, . . . , K n G {0, 1} the polynomial Kl ,... A „(O, [/) G C[C/] is of degree 
at most zero, i.e., a constant complex number, independent of U\, . . . , U n . 

Let := (^Ki,...,K n )«;i,...,K n e{o,i} an d observe that the vector 6(0, U) is a fixed point 
of the affine space A 2 ". We denote by Wl the vanishing ideal of the C-algebra C[9] 
at this point. 

Consider now the polynomial 

F := F*"' := J] (Y-(j + T Yl U^)) 

0<i<2 n -l l<i<n 

of C[T,U,Y], where denotes the i-th digit of the binary representation of the 
integer j, < j < 2 n — 1, 1 < i < n. Let q := q n . 

From the identity Y[^ l}n {Y -H{T, U, e)) = Ylo^-iiY-U+T Y[^< n U^)) 
one deduces that F is the elimination polynomial associated with the given flat 
family of zero-dimensional elimination problems G\ — 0, . . . , G n — and H. 

Let us write F q as a polynomial in the main indeterminate Y with coefficients 
tp K G C[T, U], 1 < k < 2"g, namely 

F« = F 2 ™ 9 + ^iF 2 " 9-1 + • • • + 

Observe that for 1 < k < 2 n g the polynomial (p K (0, U) G C[C/] is of degree at most 
zero. Let \ K := ip K (0, U), A := (\ K )i< K <2" q and <p := (v? K )i< K <2™< ? - Observe that A is 
a fixed point of the affine space A 2 " 9 . 

Recall that (5 is an admissible input for the procedure A and hence for that 
the final results of A^ l \(5) constitute the entries of the robust constructible map v 
defined on A ra+1 , that v represents (by means of the circuit A^ al ((3)) an admissible 
input for the procedure A^ and that 7 = As. nSL \((3) is an essentially division-free, 
parameterized arithmetic circuit with parameter domain A n+1 and input Y. 

Furthermore, recall that there exists a geometrically robust constructible map 
ip defined on the image S of v such that the entries of v* — xjj o v constitute the 
essential parameters of the circuit 7, that the entries of v and v* may be interpreted 
as polynomials of C[T, U\ and that for m being the number of components of the 
map v*, there exists a vector uj of m-variate polynomials over C such that the 
polynomials of C[T,U] which constitute the entries of oo(v*) = u o v* become the 
coefficients of the polynomial F q with respect to the main indeterminate Y. Let Tbe 
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the image of the coefficient vector 9 of H, and interpret 9 as a geometrically robust 
constructible map defined on A n+1 . Observe that T is a constructible subset of A 2 ". 
Since H is the unique final result of the circuit /3 which depends on parameters, 
we deduce from Proposition [8] that there exists a geometrically robust constructible 
map a defined on T satisfying the condition v = a o 9. This implies v* = ip o a o 9 
and, following Definition [6] (i) and |GHMSlT] . Corollary 12, that the entries of v* 
are polynomials of C[T, [/] which are integral over the local C-subalgebra C[0]gji of 
C(T,U). 

Let G C[T, [/] be such an entry. Then there exists an integer s and polynomials 
do, ai, . . . , a s G C[0] with a ^VJH such that the algebraic dependence relation 

a /i s + ai/i s_1 H + a s = (3) 

is satisfied in C[T, [/]. From ([3]) we deduce the algebraic dependence relation 

a (0, f/)/i(0, + ai(0, t/)/i(0, f/) 5 - 1 + ■ ■ ■ + a s (0, C/) = (4) 

in C[U}. 

Since the polynomials belong to C[0] and 9(0, U) is a fixed point of 

A 2 ", we conclude that «o := a (0, U), ot\ := ai(0, Z7), . . . , a s := a s (0, U) are complex 
numbers. Moreover, a VJH implies a ^ 0. 

Thus Q may be rewritten into the algebraic dependence relation 

a /i(0, U) s + ai/i(0, f/) s_1 + • • • + a s = (5) 

in C[C7] with a + 0. 

This implies that the polynomial fi(0,U) of C[Z7] is of degree at most zero. 
Therefore w := v*(0, U) is a fixed point of the affine space A m . 

Since 7 computes the polynomial F q and F q has the form F q = Y 2nq +(fiY 2nq ~ 1 + 
■ ■ ■ + if2 n q with ip K G C[T, U), 1 < k < 2 n q, we see that = (v?k)i< k <2™ (? m ay be 
decomposed as follows: 

(f = 00 (v*) = bJ O V* . 

Recall that A = (A K )i< K <2« g with A K := <p K (0, U), 1 < k < 2 n q, is a fixed point of 
the affine space A 2 ". 

For 1 < k < 2 n q we may write the polynomial (p K G C[T, U] as follows: 

V^k = A K + A K T + terms of higher degree in T (6) 

with A K G C[Z7]. From |CGH + 03] , Lemma 6 we deduce that the elimination poly- 
nomial F has the form F = Y 2 " + B{Y 2n ~ l H h -B 2 «, where for 1 < / < 2 n the 

coefficient 5; is an element of C[T,U] of the form 

B[ = (—1)' ji . . . j; + XX; + terms of higher degree in T, 

l<jl<-<jl<2 n 



45 



where L\, . . . , L 2 n G C[U] are C-linearly independent. 

Choose now different complex numbers T]i, ... , 772" from C — {j G Z; < j < 2 n } 
and observe that for 1 < k' < 2 n the identities 

— (0,U, VK ,) = qF^\0,U, VK ,) — (0,U,r ]K ,)=q J] - j)^ £ 

0<j<2™ l<K2 n 

and 

^(0,C/,„..)= A^>- 

hold. 

Since L±, . . . , L 2 n are C-linearly independent, we deduce from the non-singularity 
of the Vandermonde matrix (v^~ l )i<i,K'<2 n that 2 n many of the polynomials Ai, . . . , 
A 2 n g of C[U] are C-linearly independent. 

Consider now an arbitrary point u G A n and let e u : A 1 — > A m and 5 U : A 1 — > A 2 ™ 9 
be the polynomial maps defined for t G A 1 by e u (t) := u*(t,u) and 5 u (t) := (f(t,u). 
Then we have e„(0) = z/*(0,w) = w and S u (0) = f(0,u) = A, independently of u. 
Moreover, from ip = u o v* we deduce i u = w o e u . 

Thus ^ implies 

(A x («),..., A 2 n») = ^(0,«) = £(0) = (Du) w (e> u (0)l (7) 

where (Du) w denotes the (first) derivative of the m-variate polynomial map u at the 
point w G A m and S' U (Q) and e' u (0) are the derivatives of the parameterized curves 
5 U and e u at the point G A 1 . We rewrite now ^ in matrix form, replacing (Dtu) w 
by the corresponding transposed Jacobi matrix M G A mx2 ™ 9 and S ; u (0) and e^(0) by 
the corresponding points of A 2 " 9 and A m , respectively. 
Then ^ takes the form 

(A 1 («),...,A 2 n,(«)) = e / u (0)M, (8) 

where the complex (m x 2 ra g)-matrix M is independent of u. 

Since 2 n many of the polynomials A 1; . . . , A 2 n 6 C[[/] are C-linearly indepen- 
dent, we may choose points U\, . . . , w 2 n G A ra such that the complex (2 n x 2 n g)-matrix 

N := (A K (m)) x<i< 2 ™ 

\<K<2 n q 

has rank 2™. 

Let K be the complex {2 n x m)-matrix whose rows are e' Ul (Q), • • • , e^ 2n (0). 
Then ^ implies the matrix identity 

N — K ■ M. 
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Since N has rank 2™, the rank of the complex (m x 2 ra )-matrix M is at least 2 n . 
This implies 

m > 2 n . (9) 

Therefore the circuit 7 contains m > 2 n essential parameters. 

Let L be the number of essential multiplications executed by the parameterized 
arithmetic circuit 7 and let V be the total number of multiplications of 7, excepting 
those by scalars from C. Then, after a well-known standard rearrangement [PS73J 
of 7, we may suppose without loss of generality, that there exists a constant c > 
(independent of the input circuit 7 and the procedure A) such that L > cm 2 and 
V > cm holds. 

From the estimation ^ we deduce now that the circuit 7 performs at least 
f2(2a) essential multiplications and at least Q(2 n ) multiplications, including also 
multiplications with parameters. This finishes the proof of the theorem. ■ 



Observations Let assumptions and notations be as before. In the proof of The- 
orem 



TO] we made a substantial use of the output isoparametricity of the procedure 
A^ when we applied Proposition [8] in order to guarantee the existence of a geo- 
metrically robust constructible map a defined on T which satisfies the condition 
v = a o 9. The conclusion was that the entries of v* = ip v are polynomials of 
C[T, U] which are integral over C[#]<h. This implied finally that z/*(0,?7) is a fixed 
point of the affine space A m . Taking into account the results of |CGH + 03] , Sections 
3.2 and 5.1 it suffices to require that the procedure A^ is output coalescent in order 
to arrive to the same conclusion. This means that Theorem [10] remains valid if we 
require only that the procedure A^ is output coalescent. 



In the proof of Theorem 10 we have exhibited an infinite sequence of flat fami- 
lies of zero-dimensional elimination problems represented by robust parameterized 
arithmetic circuits of small size, such that any implementation of their associated 
elimination polynomials, obtained by a procedure of our extended computation 
model which solves the given elimination task for any instance, requires circuits 
of exponential size. 



The statement of Theorem 10 may also be interpreted in terms of a mathe- 
matically certified trade-off of quality attributes. Suppose for the moment that 
we had built our model for branching parsimonious computation in the same way 



as in Section 3.3, omitting the requirement of isoparametricity for recursive rou- 



tines, however. Recall that this requirement implies the output isoparametricity 



of any algorithm of our extended computation model and recall from Section 3.3.2 



that well behavedness under reduction is a quality attribute which implies output 



isoparametricity and therefore also the conclusion of Theorem [10 

A complexity class like "exponential time in worst case" represents also a quality 
attribute. Thus we see that the quality attribute "well behavedness under reduction" 
implies the quality attribute "exponential time in worst case" for any essentially 
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division-free procedure of our extended computation model which solves algorith- 
mically the general instance of any given flat family of zero-dimensional problems. 



The proof of Theorem 10 depends substantially on the decomposition of the 



elimination procedure A into two subprocedures A^ and A^ satisfying conditions 



(i) and (ii) of Section 3.3.3 We are now going to justify this architectural restriction 
on the procedure A for the particular case of elimination algorithms. 

As at the beginning of this section, let U := (Ui, . . . , U r ), X := (Xi, . . . ,X n ), 
Gi, . . . , G n , H e C[U, X] and F G C[U, Y] such that Gi = 0, . . . , G n = and H con- 
stitute a flat family of zero-dimensional elimination problems and F its associated 
elimination polynomial. Suppose that G\, . . . ,G n and H are implemented by an 
essentially division-free, robust parameterized arithmetic circuit /3 with parameter 
domain A r and inputs Xi, . . . , X n . 

All known algorithms which solve the general instance of any flat family of zero- 
dimensional elimination problems may be interpreted as belonging to our restricted 
set of procedures. They compute directly the elimination polynomial F (and not 
an arbitrary power of it). Thus let A be such a known algorithm and let A^- 1 ' and 
A^ be the subalgorithms which compose A in the same way as before. Then A^ 
computes the coefficients of F, where F is considered as a polynomial over C[U] 
in the indeterminate Y. The subalgorithm A^ may be interpreted as the Horner 
scheme which evaluates F from its precomputed coefficients and Y. 

Observe that F, and hence deg y F, depends only on the polynomials G\, . . . , G n 
and H, but not on the particular circuit (3. Therefore deg y F is determined by ^((3), 
where ip is the low level program of the algorithm A. 

For any parameter instance u E A r we may think (A^((3)Y U ^ as a constraint 
database (in the sense of [HK04J and [GHKa] ) which allows to evaluate the uni- 
variate polynomial F^ £ C[Y] as often as desired for arbitrary inputs y E A 1 , 
using each time a number of arithmetic operations in C, namely deg y F, which does 
not depend on the non-scalar size of j3. 

Moreover A satisfies the following condition: 

(D) There exist non-decreasing real valued functions G\ > and C% > depending 
on dynamic integer vectors, such that for Lp and being the non-scalar 

sizes of the circuits (3 andA((3), the inequality 

holds. 

Let now A be an arbitrary, essentially division-free algorithm of our extended 
computation model which solves the general instance of any flat family of zero- 
dimensional elimination problems and let /3 be an input circuit for A which rep- 
resents a particular family of such problems. Let F be the associated elimination 
polynomial. 
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Then the complexity of the algorithm A is only competitive with known elimi- 
nation algorithms if we require that the number of essential additions and multipli- 
cations of *4.fmai(/3) is bounded by 2 • deg y F. This leads us to the requirement that 
A must be decomposable in two subalgorithms A^ and A^ as above. 

Therefore any elimination algorithm of our extended computation model which is 
claimed to improve upon known algorithms for all admissible input circuits /3, must 
have this architectural structure. In particular, such an algorithm cannot call the 
input circuit /3 when the output variable Y became already involved. This justifies 



the architectural restriction we made in the statement and proof of Theorem 10 

Observe that the final results of the circuit A^'(/3) form a geometrically robust 
constructible map defined on the parameter domain of the circuit /3. For a given pa- 
rameter instance, the value of this map allows to compute the value of the coefficient 
vector of the elimination polynomial F on this instance. 

Moreover, the competitivity of A with known elimination algorithms requires 
that A must satisfy condition (D). 



From Theorem 10 and its proof we deduce now the lower bound 



m^{C 1 (^ n )),C 2 ^(P n ))} = fi(^), 

where 5 n is the geometric degree of the subvariety of A r x A ra+1 defined by the 
polynomials G^, . . . ,Gn,Y — (observe 5 n = 2 n ). Adding to (3 n a suitable 
addition node we obtain a totally division-free new circuit fl* which represents 
. . . ,Gn and Y — H^ n \ Observe that for each (s,u) G A 1 x A n the degree 
pattern of the polynomials G^, . . . , Gn, Y — H(s, u, X) is constant and the system 
degree is 5 n . The polynomial F^ is the output of the Kronecker algorithm applied to 
P* and the variable Y. Therefore the algorithm A produces on (3 n the same output 
as the Kronecker algorithm applied to 0* and the variable Y. We conclude now 
from Lp* — 0(n) that the Kronecker algorithm is nearly optimal in our extended 
computation model. 

In our computation model, algorithms are transformations of parameterized 
arithmetic circuits over one and the same parameter domain. This represents a 



substantial ingredient for the proof of Theorem 10 If we allow branchings which 



lead to subdivisions of the parameter domain of the input circuit, the conclusion of 



Theorem 10 may become uncertain (see [G HKbj ). 

Our computation model is also restrictive in another sense: 
suppose that there is given an essentially division-free, robust parameterized arith- 



metic circuit f3 n evaluating the polynomial as in the proof of Theorem 10 and 
an essentially division-free procedure B of our extended computation model which 
recomputes from the input f3 n . Then, the output of B on (3 n is an essentially 
division-free robust parameterized arithmetic circuit of size Q(2 n ), although the size 
of f3 n is 0(n). 
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4.2 The elimination of a block of existential quantifiers 



Let notations be the same as in the proof of Theorem 10 in Section |4.1[ Let n G N, 
Si, . . . , S n new indeterminates, S := (Si, . . . , S n ), G^ := X\ — Xi — Si, . . . , Gn := 
Xl - X n - S n and again ff<»> := Ei<,<„ ^X, + T Ui< t < n ^ + % - l)X t ). 

Observe that the polynomials G^, . . . , G^ form a reduced regular sequence in 
C[S, T, U, X] and that they define a subvariety V n of the affine space A n xA 1 xA"x 
A n which is isomorphic to A™ x A 1 x A n and hence irreducible and of dimension 
2n + 1. Moreover, the morphism V n -¥ A™ x A 1 x A n which associates to any point 
(s, t, u, x) G V n the point (s, t, u), is finite and generically unramified. Therefore the 
morphism 7c n : V n — > A n x A 1 x A n x A 1 which associates to any (s,t, u, x) G V n 
the point (s,t,u,H^(t,u,x)) G A n x A 1 x A n x A 1 is finite and its image n n (V n ) 
is a hypersurface of A n x A 1 x A" x A 1 with irreducible minimal equation G 
C[S, T,U,Y\. 

Hence = 0, . . . , G^ = and represent a flat family of zero-dimensional 
elimination problems whose associated elimination polynomial is just F^ n \ 
Observe that deg F^ = deg y = 2 n and that for G A n the identity 

F (n) (0, T, U, Y) = F (n) (T, U, Y) holds, 

where F^ is the elimination polynomial associated with the flat family of zero 
dimensional elimination problems given by Xf — X\ = 0, . . . , X\ — X n = and 
Since F^ is irreducible, any equation of C[S, T, U, Y] which defines 7t n (V n ) in 
A™ x A 1 x A" x A 1 is without loss of generality a power of F^ n \ 

We consider now Si, . . . , S n , T, U\, . . . , U n as basic parameters, Xi, . . . , X n as 
input and Y as output variables. 

Let A' be an essentially division-free procedure of our extended computation 
model satisfying the following condition: 

A' accepts as input any robust parameterized arithmetic circuit (3 which represents 
the general instance of a flat family of zero-dimensional elimination problems with 
associated elimination polynomial F and *4.' final (/3) has a single input Y and a single 
final result which defines the same hypersurface as F. 

With this notions and notations we have the following result. 

Proposition 11 There exist an ordinary division-free arithmetic circuit f3 n of size 
0(n) over C with inputs Si, ... , S n , T, Ui, . . . , U n , Xi, . . . , X n and final results 
G^, . . . , Gn, H( n \ The essentially division-free, robust parameterized arithmetic 
circuit % := A'fl nal ($ n ) depends on the basic parameters Si, ... , S n , T , Ui, . . . ,U n 

and the input Y and its single final result is a power of F^ . The circuit 7 n performs 
at least Q(2 2 ) essential multiplications and at least Q(2 n ) multiplications with param- 
eters. As ordinary arithmetic circuit over C with inputs Si, ... , S n , T, Ui, . . . ,U n 
and Y , the circuit j n has non-scalar size at least Q(2 n ). 

Proof. The existence of an ordinary division-free arithmetic circuit as in the state- 



ment of Proposition 11 is evident. The rest follows immediately from the proof of 
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Theorem 10 in Section 4T by restricting the parameter domain A" x A 1 x A n of /3 n 
and 7„ to A 1 x A ra , i.e., by specializing S to G A n . Observe that this restriction of 
7„ may become an inconsistent circuit, but this does not affect the argumentation 
which is based on the consideration of suitable geometrically robust constructible 
functions. ■ 

Suppose now that there is given another essentially division-free procedure A" 
of our extended computation model satisfying the following condition: 
A" accepts as input any robust arithmetic circuit /3 which represents the general 
instance of a flat family of zero-dimensional elimination problems with associated 
elimination polynomial F and there exists a Boolean circuit b in as many variables 
as the number of final results of A'l ngX ((3) such that the algebraic variety defined by 
F coincides with the constructible set which can be described by plugging into b the 
final results of A£ naA (/3) as polynomial equations. 

Observe that this represents the most general architecture we can employ to 



adapt in the spirit of Section |3.3.3| our extended computation model for functions 
to parametric decision problems. 

Let s6N and Ai,...,A s new indeterminates with A := (Ai, . . . ,A S ). We sup- 
pose that there is given an essentially division-free procedure B of our extended 
computation model which accepts as input any essentially division-free, robust pa- 
rameterized arithmetic circuit 7 with the basic parameters Ai,...,A s and the input 
variable Y, such that i3fi na i(7) represents, by its output nodes, in C[A, Y] the mul- 
tiplicative decomposition of the final results of 7 by their greatest common divisor 
and complementary factors. 

In this sense, we call the procedure B a GCD algorithm. 

Let ipA" and ipg be the given low level programs of the procedures A!' and B. We 
require that A!' and B are competitive with known algorithms which solve the same 



tasks. Following our argumentation in Section 4T we may therefore suppose that 
there exist four non-decreasing real valued functions G\ > 0, C% > and D\ > 0, 
D2 > which depend on dynamic integer vectors and which satisfy the estimates 

L A >>v» < C^i^Lp + C 2 (VU"(/?)) 

and 

£fl( 7 ) < Di(Vfe(7))£ 7 + Aj(Vte(7))- 

We consider again the ordinary division-free arithmetic circuit /?„ of Proposition 
^ which represents the polynomials G^\ . . . , G^ and H^ n \ 

With these notions and notations we may now formulate the following statement. 

Theorem 12 Let assumptions and notations be as before. Then we have 

,2i 



max{C i ^ A 4P n )),D t (M^ nal (Pn)))^ = 1,2} = Q(- 



n 
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Proof. If we plug into the Boolean circuit b the final results of A'^ a&l {Pn) as poly- 
nomial equations, we obtain by assumption a description of the hypersurface n(V n ) 
of the affine space A n x A 1 x A" x A 1 . This implies that between the final results 
of ^4.fi na i(/3n) there exists a selection, say the polynomials Pi, . . . , P m and Pi, . . . , R t 
of C[S, T, U, Y] such that the formula 

P 1 = 0A---AP m = 0Ai?i^0A---A J R 4 ^0 

defines a nonempty Zariski open (and dense) subset of the irreducible surface fr(V n ) 
of A n x A 1 x A n x A 1 . 

Let R := R\ . . . R t and observe that the greatest common divisor of Pi, ... , P m 
has the form (p( n )) 9 ■ Q, where q belongs to N and Q is the greatest common divisor 
of Pi, ... , P m , R. Therefore we may compute (p( n ))<? in the following way: erasing 
suitable nodes from the circuit A'^ nal (f3 n ) and adding t — 1 multiplication nodes 
we obtain two robust parameterized arithmetic circuits 7i™' ) and 73 with basic 
parameters Si, ... , S n , T, Ui, . . . , U n and input Y whose final results are Pi, ... , P m 
and Pi, ... , P m , R respectively. 

Between the final results of i3 nna i(7j™' ) ) and Bf[ nSi \('y^) are the polynomials 
(F^) q -Q and Q. Applying the procedure B to the union of -Bfinai^i^) and -Bfinai^"^) 
we obtain finally an essentially division-free, robust parameterized arithmetic circuit 
with basic parameters Si, . . . , S n , T, U\, . . . , U n and input Y whose single final result 
is (F n ) q . 

Joining the circuits A"((3 n ), i3 nna i(7i ?1 ' > ), Banaiili^) an d the final division node 
we obtain an ordinary arithmetic circuit of non-scalar size at most 

1 + 3L B(A>' nal n )) ^ 

1 + 3(D 1 ^ B (A^ nal n )))L A , n&i n) + DMA'LM))) < 

3C 2 (^(A„))Pl(^Knal(/3n))) + D 2 {fa (^al A) ) ) ■ 

On the other hand we deduce from Theorem | 

L- pn =0{n) and 1+3^^^ = 0(2"). 



This implies the estimate of Theorem 12 



In a simple minded understanding, Theorem 12 says that in our extended com- 
putation model either the elimination of a single existential quantifier block in a 
prenex first-order formula of the elementary language of C or the computation of 
a greatest common divisor of a finite set of circuit represented polynomials requires 
exponential time. Complexity results in this spirit were already obtained in [GTJOlj 



and jCGH + 03] (compare also Proposition 11 and Observation in Section 4.1). 
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4.3 Arithmetization techniques for Boolean circuits 

Let m G N and let 0,1 and Z\,...,Z m be given constants and variables. Let 
Z := {Zxi • • • , Z m ). Following the context we shall interpret 0, 1 as Boolean values 
or the corresponding complex numbers and Z\, . . . , Z m as Boolean variables or inde- 
terminates over C. With A, V, ~ we denote the Boolean operations "and", "or" and 
"not". A Boolean circuit b with inputs Z%, . . . , Z m is a DAG whose indegree zero 
nodes are labelled by 0, 1 and Z 1; . . . , Z m and whose inner nodes have indegree two 
or one. In the first case an inner node is labelled by A or V and in the second by ~. 
Some inner nodes of b become labelled as outputs. We associate with b a semantics 
as follows: 

- indegree zero nodes which are labelled by 0, 1 become interpreted by the cor- 
responding constant functions {0, l} m — » {0, 1}, 

- indegree zero nodes which are labelled by Z%, . . . , Z m become interpreted by 
the corresponding projection function {0, l} m — > {0, 1}, 

- let p be an inner node of b of indegree two whose parent nodes p\ and p2 are 
already interpreted by Boolean functions g Pl ,g P2 '■ {0, l} m —> {0,1}. If p is 
labelled by A, we interpret p by the Boolean function g p := g pi A g P2 and if p 
is labelled by V, we interpret p by the Boolean function g p := g pi V g P2 , 

- let p be an inner node of b of indegree one whose parent node p' became already 
interpreted by a Boolean function g p i : {0, l} m — > {0, 1}. Then we interpret p 
by the Boolean function g p := g~^. 

For a node p of b we call g p the intermediate result of b at p. If p is an output 
node, we call g p a final result of b. 

An arithmetization /3 of the Boolean circuit b consists of the same DAG as b 
with a different labelling as follows. 

Let U, V be new indeterminates over C. The constants 0, 1 become interpreted 
by the correspondent complex numbers and Z\, . . . ,Z m as indeterminates over C. 
Let p be an inner node of (3. If p has indegree two, then p becomes labelled by a 
polynomial R p G Z[U, V] and if p has indegree one by a polynomial R p G Z[Z7]. The 
output nodes of and b are the same. 

Representing for each inner node p of /3 the polynomial G p by a division-free 
ordinary arithmetic circuit over Z in the inputs U, V or U, we obtain an ordinary 
division-free arithmetic circuit over Z in the inputs Zi, . . . , Z m . 



Just as we did in Section 3.3.2 we may associate with (3 a semantics which 
determines for each node p of (5 a polynomial G p G Z[Z]. We say that /3 is an 
arithmetization of the Boolean circuit b if the following condition is satisfied: 
for any node p of b and any argument z G {0, l} m the Boolean value g p (z) coincides 
with the arithmetic value G p {z) (in a somewhat imprecise notation: g p (z) = G p {z)). 

An example of an arithmetization procedure is given by the map which associates 
to each node p of b a polynomial [g p ] of 7L\Z\ satisfying the following conditions: 
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[0] = 0, [1] = 1, \Z X \ = Z U ..., [Z m ] = z 



rn 



- for p an inner node of indegree two of b with parents p\ and p2- 



[9 P ] = [9 Pl ] 



[g P2 ] if the label of p is A 



and 



l9 P ] = [g Pl ] + [g P2 ] 



[dpi] ' \9p2\ if the label of p is V 



- for p an inner node of indegree one of b with parent p'\ 



[9p\ = 1 - lap'}- 



The resulting arithmetic circuit is called the standard arithmetization of b (see, 
e.g., |Sha92j and [EF91 J). 

Let n, r G N and U\, . . . , U r , Xi, . . . , X n be new variables. For m := n + r we 
replace now Z by U and X, where U := (E/i, . . . , U r ) and X := (X 1; . . . , X n ). We 
shall interpret U\, . . . ,U r as parameters and Xi, . . . , X n as input variables. 

Let b be a Boolean circuit with the inputs Ui, . . . , U r , Xi, . . . , X n and just one 
final result h : {0, l} r x {0, 1}" {0, 1}. 

We wish to describe the set of instances u G {0, l} r where h(u, X 1 , . . . , X n ) is a 
satisfiable Boolean function. 

For this purpose let us choose an arithmetization (3 of b. We interpret f3 as an 
ordinary arithmetic circuit over Z with the parameters U\, . . . ,U r and the inputs 
Xi, . . . , X n . The single final result of /3 is a polynomial G Z[Z7, X] which satisfies 
for any u G {0, l} r , x G {0, l} n the following condition: 



Without loss of generality we may suppose that the polynomials X\ — X 1; . . . , 
X^ — X n are intermediate results of /3. We relabel now (3 such that these polynomials 
and H become the final results of (3. Observe that X\ — X\ = 0, . . . , X\ — X n = 
and H represent a flat family of zero-dimensional elimination problems. 

Let Y be a new indeterminate and let F G Z[£7, Y] the associated elimination 
polynomial. One verifies easily the identity 



Let A be an essentially division-free procedure of our extended computation 
model which solves algorithmically the general instance of any flat family of zero- 
dimensional elimination problems. Then j3 is an admissible input for A and there 
exists an integer q G N such that F q is the final result of A& ng i(f3). 



We consider now the task to count for any u G {0, l} r the number k of instances 
x G {0, l} n with h(u,x) = 1. 



h(u, x) = H(u, x). 



F(U,Y)= ][ (Y-H(U,x)). 



xe{o,i}" 



54 



The polynomial F q encodes two possible solutions of this task. 

The first solution is the following: let I be the order of the univariate polynomial 
F q (u, Y) at zero. Then q divides I and we have k = 2 n — -. 

The second and more interesting solution is the following: write F q = Y 2nq + 
(p 1 Y 2Uq ~ 1 + ■ ■ • + if2^ q with <fi, . . . , ip 2 n q G Z[U]. Then ifi(u) is an integer which is 
divisible by q and we have k = — . 

Observe also deg ipi < deg v H. 

These considerations show the relevance of an efficient evaluation of the poly- 
nomial F q (e.g., by the circuit A& na .i(f3)) ■ 

We ask therefore whether A^ na \(j3) can be polynomial in the size of the Boolean 
circuit b. The following example illustrates that the answer may become negative. 

In the sequel we are going to exhibit for r := 2n + 1 a Boolean circuit b of size 
0{n) which evaluates a function h : {0, l} r x {0, l} n — > {0, 1} such that the stan- 
dard arithmetization (3 of b represents a flat family of zero-dimensional elimination 
problems with associated elimination polynomial F and such that any essentially 
division-free procedure A of our extended computation model that accepts the in- 
put (3 and computes by means of A& na i(f3) a power of F, requires time fl(2 n ) for this 
task. This means that it is unlikely that algorithms designed following the paradigm 
of object-oriented programming are able to evaluate the polynomial <pi efficiently 

On the other hand, since the degree of <pi is bounded by deg^ H and therefore 
"small", there exists a polynomial time interactive protocol which checks for any 
u G {0, l} r and any c G Z the equation ifi(u) = c. Thus this problem belongs to 
the complexity class IP (see [LFKN92J for details). 

We are now going to exhibit an example of a Boolean circuit which highlights 
the unfeasibility of our computation task. 

For this purpose let r := 2n + 1 and Si, . . . , S n , T, Ui, . . . , U n parameters and 
Xi, . . . , X n input variables and let S := (Si, . . . , S n ) and U := (Ui, . . . , U n ). 

We consider the Boolean function h : {0, l} 2n+1 x {0, l} n — > {0, 1} defined by 
the Boolean formula 

0:= /\ (JQV0%AJQ))V(TA /\ (XlviUiAXi))). 

l<i<n l<i<n 

From the structure of the formula we infer that h can be evaluated by a Boolean 
circuit b of size 0{n) in the inputs Si, . . . , S n , T, U\, . . . , U n . 

Let j3 be the standard arithmetization of the Boolean circuit b and let H be the 
final result of (3. Observe that the total, and hence the non-scalar size of (3 is 0(n). 
Then we have 

h= n w + (i - n (i + (Si - i)Xi))T n a + (^ - i^. 

l<i<n l<i<n l<i<n 

Observe that the equations X\ — X\ = 0, . . . , X\ — X n = and the polynomial 
H represent a flat family of zero-dimensional elimination problems. Let F be the 
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associated elimination polynomial. Then F can be written as 



F = Y 2n +B l Y 2n - 1 + ...+B 2n = J] (Y-( J] S^+(l- J] S^)T n U^)) 

0<j<2 n l<i<n l<i<" l<i<n 

with 

B k = (-i) fc e n ( n *p + (i - n n ^ 

0<ji<---<jfc<2" l</i<fc l<i<n l<i<n \<i<n 

for 1 < k < 2 n . 
Let 

^ ■= e En 5 P c 1 - n n ^ n ^ 

0<ji< -<j fc <2" l<fo<fe l<i<n l<i<n l<i<n l<i<n 

where 1 < k < 2 n . 
Then we have 

B k = (-l) k E II ^ ' ' ' II S i kU + L k-T+ terms of higher degree in T 

0<ji< --<j fc <2" l<i<n l<i<n 

Let e : A 2 ™ — > A 2 ™ be the morphism of affine spaces which assigns to each point 
z G A 2 ™ the values of the elementary symmetric functions in 2 n variables at z. 
Observe that the Jacobian of e at (rii<j<n >S'' J ' i )o<i<2™ is a non-singular (2™ x 2 n )- 
matrix N(S). The polynomials 1 < k < 2 n are obtained by applying N(S) to 

((l-n«i<„ Ui< t < n U? ]i ) < j<2n . Since the monomials Ui<i<n U?\ 0<j< 2", 
are linearly independent over C(S) we conclude that the polynomials L^, 1 < k < 2 n 
have the same property. 

With this preparation we are now able to repeat textually the arguments in the 
proof of Theorem 10 of Section 4^ in order to show the following statement. 

Theorem 13 Let assumptions and notations be as before and let A be an essen- 
tially division free procedure of our extended computation model which accepts the 
arithmetic circuit (3 as input. Suppose that Afi na i((3) has a unique final result and 
that it is a power of the elimination polynomial F. Then the non-scalar size of 
Afi na i{/3) is at least n(2 n ). 



4.4 The multivariate resultant 

Let Ui, . . . , U m be basic parameters and let X 1: . . . , X n or X , . . . , X n be input vari- 
ables. Let Gi, . . . , G n , H G C[U, X] be a flat family of zero-dimensional elimination 
problems such that for any u G A m the homogenizations of Gi(u, X), . . . , G n (u, X) 
(by Xo) have no common zero at infinity. Let F G C[U, Y] be the corresponding 
elimination polynomial and let R G C[U] be the (multivariate) resultant of the 
homogenizations of G\, . . . , G n , H. Then we have R = F(U, 0). 
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On the other hand, F is the resultant of the homogenizations of Gi, . . . ,G n 
and H — Y. Thus, multihomogeneous resultants and elimination polynomials of 
flat families of zero-dimensional elimination problems are closely related from the 
algebraic point of view. 

From the computational point of view this relation is more intricate. 

To the degree pattern of the homogenizations of Gi, . . . , G n , H there corresponds 
a generic resultant. We may take a computation of this resultant and specialise its 
inputs to the coefficients of Gi, . . . , G n , H — Y with respect to Xi, . . . , X n . If this 
specialized computation can be simplified by means of reductions, we may expect 
to gain something. Proceeding in this way we obtain an algorithm which may 
be interpreted as an elementary routine of our computation model. Observe that 
this elementary routine does not use joins of two subroutines that contain each a 
recursion. 

We are going to show that reductions do not produce a general improvement of 
traditional resultant computations in the sense described above. 

Theorem 14 Consider the following, with respect to Xq, . . . , X n homogeneous, poly- 
nomials 

X\ - X Xx, ...,Xl- X X n , YX£ - J2 ^X^Xi - T J] (X + ([/, - 1)X,), 

l<i<n l<i<n 

which are supposed to be given by an essentially division-free arithmetic circuit (3 in 
the basic parameters T, U\, . . . , U n , Y and the input variables Xq, . . . , X n . Suppose 
furthermore that Y is the last basic parameter introduced by (3 and observe that such 
a circuit (3 of size 0(n) exists. Let A be an essentially division-free elementary 
routine of our computation model which on input (3 evaluates the resultant of the 
polynomials above with respect to the variables Xq, . . . , X n . Suppose that A does not 
use joins of two subroutines that contain each a recursion. Then the output circuit 
Afinai((3) has size at least f2(2 n ). 

Proof. Taking into account that Y is the last basic parameter introduced by 
(3, that all recursive subroutines of A are isoparametric and that A does not use 
joins of two subroutines that contain each a recursion, we may decompose A in two 
subroutines A^ and A^ such that A^ admits (3 as input and produces a vector 
of parameters in T, Ui, . . . , U n which themselves constitute the inputs of A^ (the 
argument is somewhat tedious and is not given here). Only the routine A^ intro- 
duces the parameter Y. Therefore A = (A^^A^) constitutes a procedure of our 
computation model with basic parameters T, U\, . . . , U n , input variables X , . . . , X n 
and output variable Y . Observing that the resultant of the given polynomials equals 
the elimination polynomial of the flat family of zero-dimensional elimination prob- 
lems given by X\ - X x , . . . , X 2 n - X n , £ 1<i<n T~ l X t + T fli^ 1 + ( U i - we 



conclude from Theorem 10 that the size of Abj&x(/3) is at least f2(2 n ). 
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4.5 A family of hard elimination polynomials 

As a major result of this paper, we are now going to exhibit an infinite family of pa- 
rameter dependent elimination polynomials which require essentially division-free, 
robust parameterized arithmetic circuits of exponential size for their evaluation, 
whereas the circuit size of the corresponding input problems grows only polynomi- 
ally. This result is valid without any architectural assumption on the algorithm that 
computes these elimination polynomials. 

Let notations be as before and consider again for given n G N the polynomial 



H(n) ■= Ei< i < n 2 ^ 1 ^ + r ni< i <„(l + (^-l)^) of SectionO Observe that 
can be evaluated using n — 1 non-scalar multiplications involving X\, . . . , X n . 

The set O := {Ei<><„ + t Ui<i< n ^ + fa - "i, ■ ■ ■ , O e A« +1 } 

is contained in a finite-dimensional C-linear subspace of C[X] and therefore O and 
its closure O are constructible sets. 

From |GHMSlI) . Section 3.3.3 we deduce the following facts: 
there exist K := lQn 2 + 2 integer points £i, G Z n of bit length at most 4n 

such that for any two polynomials f,g & O the equalities /(£&) = <?(6fc)> 1 < A; < if, 
imply f = g. Thus the polynomial map H : (9 — > A x defined for / G O by 
H(/) := (/(£i), * s injective. Moreover .M := S(O) is an irreducible 

constructible subset of A^ and we have At = Finally, the constructible map 

:= H -1 , which maps .M onto (9 and .M onto O, is a restriction of a geometrically 
robust map and therefore by Corollary [6] itself geometrically robust. 

For e G {0, l} n we denote by e the map .M — > A 1 which assigns to each point v G 
jM the value 0(f) (e). From Corollary [6] we conclude that e is a geometrically robust 
constructible function which belongs to the function field C(Ai) of the irreducible 
algebraic variety Ai. 

Observe that for t G A 1 and u G A n the identities (p t (E(H^(t, u, X))) = 
(j)(E(H^(t,u,X)))(e) = ((S" 1 oE)(H^(t,u,X)))(e) = H^ n \t,u,e) hold. 

Let : = riee{o i} n ( Y ~ Then is a geometrically robust constructible 
function which maps Ai x A 1 (and hence x A 1 ) into A 1 . Considering again the 
elimination polynomial = U 0<j <2^-i( Y ~ 0' + ^Iii<*<. U h) = U^oM Y ~ 



HW(T, U, e)) of Section 4.1 we have for t G A 1 and u G A n the identities 



F^(s(if(")(t, M ,x)),r)= J] (y-0 e (s(HW(t,tt,x)))) = 

ee{0,1} " no) 

J] (Y-H^(t,u,e)) = F^(t,u,Y) 

e£{0,l} n 

Let Si, ... , Sk be new indeterminates and observe that the existential first order 
formula of the elementary theory of C, namely 

(3X 1 ) . . . (3X n )(3T)(3f/ 1 ) . . . {3U n ){Xl - X 1 = A ■ ■ ■ A X 2 n - X n = OA 
/\ Sj = # (n) (T, U,Q A Y = H {n \T, U,X)) (11) 

l<j<K 
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describes the constructible subset {(s,y) E A K+1 ;s E M.,y E A 1 , F^ n \s, y) = 
0} of A K+1 . Moreover, F^ is the greatest common divisor in C(.M)[F] of all 
polynomials of C[.M] \Y] which vanish identically on the constructible subset of A^ +1 



defined by the formula ( 11 ). Hence E C(JA)\Y] is a (parameterized) elimination 
polynomial. 

Observe that the polynomials contained in the formula ( |TT) can be represented 
by a totally division-free arithmetic circuit /3 n of size 0(n 3 ). Therefore, the formula 



(11) is also of size 0{n 



Theorem 15 Let notations and assumptions be as before and let 7 be an essentially 
division- free, robust parameterized arithmetic circuit with domain of definition M. 
such that 7 evaluates the elimination polynomial F^ n \ 

Then 7 performs at least f2(2a ) essential multiplications and at least Q(2 n ) mul- 
tiplications with parameters. 

Proof. Let 7 be as in the statement of the theorem. Without loss of generality 
we may assume that 7 has a single output node which evaluates the polynomial 
_p( n ). There exists a totally division-free arithmetic circuit 7„ of size 0(n 3 ) which 
computes at its output nodes the polynomials H^ n \T, U, 1 < k < K. 



From (10) we deduce that the join 7 * 7„ of the circuit j n with the circuit 7 at 
the basic parameter nodes of 7 is an essentially division-free, robust parameterized 
arithmetic circuit which evaluates the elimination polynomial F^ . Observe that the 
outputs of 7„ are only parameters and that only the circuit 7 introduces the variable 
Y. Moreover, there exists an isoparametry between H^ n > and the outputs of %. We 
may therefore think that the circuit 7 * j n is produced by an essentially division-free 
procedure of our extended computation model which becomes applied to the circuit 



(3 n . From Theorem 10 and its proof we deduce now that 7 * 7 n contains at least 
Q(2 n ) essential multiplications and at least Q(2 n ) multiplications with parameters. 
Since the size of 7„ is 0(n 3 ), we draw the same conclusion for 7. ■ 



Theorem 15 is essentially contained in the arguments of the proof of [GH01J, 
Theorem 5 and |CGH+03j . Theorem 4. 

Observe that a quantifier-free description of M. by means of circuit represented 
polynomials, together with an essentially division-free, robust parameterized arith- 
metic circuit 7 with domain of definition Ai, which evaluates the elimination poly- 
nomial F^ n \ captures the intuitive meaning of an algorithmic solution of the elimina- 
tion problem described by formula (11), when we restrict our attention to solutions 
of this kind and minimize the number of equations and branchings. In particular, 
the circuit 7 can be evaluated for any input point (s,y) with s E M. and y E A 1 
and the intermediate results of 7 are polynomials of C(.M)[Y] whose coefficients are 
geometrically robust constructible functions defined on M.. 
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4.6 Divisions and blow ups 



We are now going to analyze the main argument of the proof of Theorem 10 from a 
geometric point of view. 

We recall first some notations and assumptions we made during this proof. 

With respect to the in determinates X\, . . . ,X n , we considered the vector 9 of 
coefficients of the expression 



l<i<n l<i<n 

as a polynomial map A n+1 — > A 2 ™ with image T . Recall that T is an irreducible 
constructible subset of A 2 ". 

Further, with respect to the indeterminate Y, we considered the vector (p of 
nontrivial coefficients of the monic polynomial 

f= n ( y -u+ t n u ^ 

I<i<2™-1 l<i<n 

also polynomial map A" +1 — > A 2 ™. 

One sees immediately that there exists a unique polynomial map rj : T — > A 2 " 
such that if = rjoO holds. Using particular properties of 9 and <p we showed implicitly 



in the proof of Theorem 10 that rj satisfies the following condition: 



Let m be a natural number, ( : T — > A m a geometrically robust con- 
structible and 7r : A m — > A 2 ™ a polynomial map such that 77 = tc o ( holds. 
Then the condition 

m > 2" 

is satisfied. 

This means that the following computational task cannot be solved efficiently: 

Allowing certain restricted divisions, reduce the datum 9 consisting of 2 n entries 
to a datum £ consisting of only m < 2 n entries such that the vector rj still may 
be recovered from ( without using any division, i.e., by an ordinary division-free 
arithmetic circuit over C. 

Here the allowed divisions involve only quotients which are geometrically robust 
functions defined on T. 

In order to simplify the following discussion we shall assume without loss of 
generality that all our constructible maps have geometrically robust extensions to 

r. 

Let / and g be two elements of the coordinate ring C[7~] of the affine variety T 
and suppose that g 7^ holds and that the element ^ of the rational function field 

C(T) may be extended to a robust constructible function defined on T, which we 
denote also by ^, since this extension is unique. 
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Then we have two cases: the coordinate function g divides / in C[7~] or not. In 
the first case we may compute ^, by means of an ordinary division-free arithmetic 

circuit over C, from the restrictions to T of the canonical projections A 2 " — > A 1 . 
Thus ^ belongs to the coordinate ring C[T]. In the second case this is not any more 

true and C[7~][^] is a proper extension of C[7~] in C(T). In both cases C[7~][£] is 
the coordinate ring of an affine chart of the blow up of C[7~] at the ideal generated 
by / and g. We refer to this situation as a division blow up which we call essential 
if - does not belong to C[7~]. 



Therefore we have shown in the proof of Theorem [TO] that essential division blow 
ups do not help to solve efficiently the reduction task formulated before. 

A similar situation arises in multivariate polynomial interpolation (see [GHMSllJ, 
Theorem 23). 

Following |Har92] . Theorem 7.2.1 any rational map may be decomposed into a 
finite sequence of successive blow ups followed by a regular morphism of algebraic 
varieties. Our method indicates the interest to find lower bounds for the number of 
blow ups (and their embedding dimensions) necessary for an effective variant of this 
result. 

Problem adapted methods for proving lower bounds for the number of blow 
ups necessary to resolve singularities would also give indications which order of 
complexity can be expected for efficient desingularization algorithms (see [EVOOj ). 
At this moment only upper bound estimations are known |Bla09| . 

4.7 Comments on complexity models for geometric elimi- 
nation 

4.7.1 Relation to other complexity models 

The question, whether P ^ NP or Pq ^ NP^ holds in the classical or the BSS 
Turing Machine setting, concerns only computational decision problems. These, on 
their turn, are closely related to the task to construct efficiently, for a given prenex 
existential formula, an equivalent, quantifier free one (compare [BSS89] . [HM93J, 
[SS95] and [BCSS98]). In the sequel we shall refer to this and to similar, geometri- 
cally motivated computational tasks as "effective elimination" . 



Theorem 10 in Section 4.1 does not establish a fact concerning decision problems 
like the Pc ^ NPc question. It deals with the evaluation of a function which assigns 
to suitable prenex existential formulas over C canonical, equivalent and quantifier- 
free formulas of the same elementary language. 

Theorem [10] says that in our computation model this function cannot be evalu- 
ated efficiently. If we admit also non-canonical quantifier-free formulas as function 
values (i.e., as outputs of our algorithms), then this conclusion remains true, pro- 
vided that the calculation of parameterized greatest common divisors is feasible and 
efficient in our model (see [CGH+03j . Section 5). 

It is not clear what this implies for the Pc 7^ NPc question. 
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Intuitively speaking, our exponential lower complexity bound for effective geo- 
metric elimination is only meaningful and true for computer programs designed in 
a professional way by software engineers. Hacker programs are excluded from our 
considerations. 

This constitutes an enormous difference between our approach and that of Turing 
machine based complexity models, which always include the hacker aspect. There- 
fore the proof of a striking lower bound for effective elimination becomes difficult in 
these models. 

Our argumentation is based on the requirement of output parametricity which 
on its turn is the consequence of two other requirements, a functional and a non- 
functional one, that we may employ alternatively. More explicitly, we require that 
algorithms (and their specifications) are described by branching parsimonious as- 
serted programs or, alternatively, that they behave well under reductions (see Sec- 
tions 



3.3.2 and 3.3.3). 



Let us observe that the complexity statement of Theorem 10 refers to the rela- 
tionship between input and output and not to a particular property of the output 
alone. In particular, Theorem 10 does not imply that certain polynomials, discussed 
below, like the permanent or the Pochhammer polynomials, are hard to evaluate. 



Let notations and assumptions be as in Section |4.1| There we considered for 
arbitrary neN the flat family of zero dimensional elimination problems 



given by 
and 



= 0, . . . , = 0, # (n) 

z^M ._ v 2 y r<(. n ) ■— Y 2 Y 

H [n) 2 i ~ x X i + T H (1 + (Ui - l)Xt 

l<i<n Ki<n 



Let X n+ i, . . . ,Xz n _i be new indeterminates and let us consider the following poly- 
nomials 

G<gi := Xn+i -2X 2 -X U ..., Gf := X 3 - A>_i - 2^ n X^ n+u n+2<j<2n-l, 



(*2n '■— ~ U\X\ + X 1 — 1, 



G^ :— X k — Uk-2n+lXk~lX k ^2n+l + X k _iX k _ 2 n+l — -Xjfc-lj 2n + 1 < k < 3n — 1 

and 

£(n) ._ X 2n _i + TX^n^i. 

One verifies easily that G± = 0, . . . , G^'_j = 0, is another flat family of zero 
dimensional elimination problems and that both families have the same associated 
elimination polynomial, namely 

\<i<n 



F^:=H(Y-(j+T ][ U^)) 
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Suppose now that there is given an essential division-free procedure A of our ex- 
tended computation model which solves algorithmically the general instance of any 
given flat family of zero-dimensional elimination problems. 

Let n and /5* be two essentially division-free, robust parameterized arithmetic 
circuits which implement the first and the second flat family of zero dimensional 
elimination problems we are considering. 

Then (3 n and /3* are necessarily distinct circuits. Therefore *4fmai(/3n) and Aa na \(f3^) 
represent two implementations of the elimination polynomial by essentially 
division-free, robust parameterized arithmetic circuits. 

From Theorem 10 and its proof we are only able to deduce that the circuit 
-4finai(/3n) has non-scalar size at least fi(2 n ), but we know nothing about the non- 
scalar Size Of -4final(/^). 

In the past, many attempts to show the non-polynomial character of the elimina- 
tion of just one existential quantifier block in the arithmetic circuit based elementary 
language over C, employed the reduction to the claim that an appropriate candidate 
family of specific polynomials was hard to evaluate (this approach was introduced 
in |HM93j and became adapted to the BSS model in |SS95j ). 

The Pochhammer polynomials and the generic permanents discussed below form 
such candidate families. 

In Section 4.5 we exhibited a certified infinite family of parameter dependent 
elimination polynomials which require essentially division-free, robust parameterized 
arithmetic circuits of exponential size for their evaluation, whereas the circuit size 
of the corresponding input problems grows polynomially. 

Here the requirement of robustness modelizes the intuitive meaning of an algo- 
rithmic solution with few equations and branchings of the underlying elimination 
problem. 



4.7.2 The hacker aspect 

Let us finish this section with a word about hacking and interactive (zero-knowledge) 
proofs. 

Hackers work in an ad hoc manner and quality attributes are irrelevant for 
them. We may simulate a hacker and his environment by an interactive proof system 
where the prover, identified with the hacker, communicates with the verifier, i.e., the 
user of the hacker's program. Thus, in our view, a hacker disposes over unlimited 
computational power, but his behaviour is deterministic. Only his communication 
with the user underlies some quantitative restrictions: communication channels are 
tight. Hacker and user become linked by a protocol of zero-knowledge type which 
we are going to explain now. 

Inputs are natural numbers in unary representation. Each natural number rep- 
resents a mathematical object belonging to a previously fixed abstract data type of 
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polynomials. For example iieN may represent the 2 n -th Pochhammer polynomial 

?-■■= n ( T ~ti 

0<j<2 n 

or the n-th generic permanent 

Perm n := ^ ^ir(i), • • • , ^(n): 

T6Sym(n) 

where T and X n , . . . , X nn are new indeterminates and Sym(n) denotes the symmet- 
ric group operating on n elements. 

On input n e N the hacker sends to the user a division-free labelled directed 
acyclic graph r n (i.e., a division-free ordinary arithmetic circuit over Z) of size 
and claims that T n evaluates the polynomial represented by n. 

The task of the user is now to check this claim in uniform, bounded probabilistic 
or non-uniform polynomial time, namely in time n°^\ 

In the case of the Pochhammer polynomial and the permanent a suitable protocol 
exists. This can be formulated as follows. 

Proposition 16 The languages 

£poch ■= {{n, (Tj)o<j<n); n G N, 1^- is for < j < n 

a division-free labelled directed acyclic graph evaluating T— } 

and 

£perm '■= {{ n ^)'> n ^ N ,T is a labelled directed acyclic graph evaluating Perm n } 

belong to the complexity class BPP and hence to P/poly (here n e N is given in 
unary representation). 

Proof. We show only that £p oc h belongs to the complexity class P/poly. The 
proof that £p oc h belongs to BPP follows the same kind of argumentation and will 
be omitted here. The case of the language £p er m can be treated analogously and we 
shall not do it here (compare |KI04] . Section 3). 

Let n EN and let T be a division-free labelled directed acyclic graph with input 
T and a single output node. Let T' be the division-free labelled directed acyclic 
graph which is given by the following construction: 

- choose a labelled acyclic graph fi n of size n + 0(1) with input T and with 
T — 2 2 " as single final result 

- take the union Y of the circuits Y and r * \i n and connect the two output nodes 
of T by a multiplication node which becomes then the single output node of 
the resulting circuit Y' . 
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From the polynomial identity T— = T?—^(T) • T^(T-2 2 " ) one deduces easily 
that r" computes the polynomial T— if and only if T computes the polynomial T— — . 

For < j < n let Tj be a division-free labelled directed acyclic graph with input 
T and a single output node. 

Suppose that in the previous construction the circuit T is realized by the labelled 
directed acyclic graph T n _i. Then one sees easily that (n, (Tj)o<j< n ) belongs to £p oc h 
if and only if the following conditions are satisfied: 

(i) the circuit T computes the polynomial T, 
(ii) the circuits T' and T n compute the same polynomial, 
(Hi) (n - 1, {Tj)o<j< n -i) belongs to £ Poc h- 

Therefore, if condition (ii) can be checked in non-uniform polynomial time, the 
claimed statement follows. 

For < j < n let Lj and L be the sizes of the labelled directed acyclic graphs 
Tj and T' and observe that L = 2L„_i + n + 0(1) holds. 

Let P n -i and P be the final results of the circuits T n _i and V . From [CGH + 03| . 
Corollary 2 we deduce that there exist m := 4(L + 2) 2 + 2 integers 71, . . . , 7 m e Z 
of bit length at most 2(L + 1) such that the condition (ii) above is satisfied if and 
only if 

(iv) P„-i( 7 l) = P(7l), ■ ■ ■ , Pn-l{l m ) = P{lm) 

holds. 

From [HMJ we infer that condition (iv) can be checked by a nondeterministic 
Turing machine with advise in (non-uniform) time 0(L 3 ) = 0((L n _i + n) 3 ). 

Applying this argument recursively, we conclude that membership of (n, (^j)o<j<n) 
to £poch may be decided in non-uniform time 0(J2o<j< n (Lj + j) 3 ) and therefore in 
polynomial time in the input size. Hence the language £p oc h belongs to the complex- 
ity class P/poly. The proof of the stronger result, namely £p oc h £ BPP, is similar. 
■ 

Finally we observe that for n G N the Pochhammer polynomial T 2 ™ is the asso- 
ciated elimination polynomial of the particular problem instance^ given by T := 0, 
of the flat family of zero-dimensional elim ination problems = 0, . . . , = 
0,H^ n \ which we considered in Section 4.1 

From the point of view of effective elimination, the sequence of Pochhammer 
polynomials becomes discussed in [HM93J (see also |SS95j ). From the point of view 
of factoring integers, Pochhammer polynomials are treated in |Lip94| . 

4.7.3 Final comment 

Let us mention that our approach to effective elimination theory, which led to The- 
orem 



10 and preliminary forms of it, was introduced in [HMPW98J and extended in 



[GH01J and |CGH+03j . 
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The final outcome of our considerations in Sections 4A and 4/7 is the following: 
neither mathematicians nor software engineers, nor a combination of them will ever 
produce a practically satisfactory, generalistic software for elimination tasks in Alge- 
braic Geometry. This is a job for hackers which may find for particular elimination 
problems specific efficient solutions. 
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